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Keldysh representation of the functional integral for the interacting electron system with disorder 
is used to derive microscopically an effective action for dirty superconductors. In the most general 
case this action is a functional of the 8x8 matrix Q(t, t') which depends on two time variables, and on 
the fluctuating order parameter field and electric potential. We show that this approach reproduces, 
without the use of the replica trick, the well-known result for the Coulomb-induced renormalization 
of the electron-electron coupling constant in the Cooper channel. Turning to the new results, we 
calculate the effects of the Coulomb interaction upon: i) the subgap Andreev conductance between 
superconductor and 2D dirty normal metal, and ii) the Josephson proximity coupling between 
superconductive islands via such a metal. These quantities are shown to be strongly suppressed 
by the Coulomb interaction at sufficiently low temperatures due to both zero-bias anomaly in the 
density of states and disorder-enhanced repulsion in the Cooper channel. 
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I. INTRODUCTION 

Electron transport in hybrid superconductive-normal systems at low temperatures is governed by the Andreev 
reflection. Both finite- voltage conductance G a between superconductive and normal electrodes, and Joscphson critical 
current I c between two superconductive banks, separated by normal region, are determined by the Cooper pair 
propagation in the normal metal. The theory of the Andreev conductance (without Coulomb effects) was developed 
in, e. g., whereas the Josephson coupling was calculated (with the account of short-range electron interaction) 

in Q . When normal conducting region is made of a dirty metal film, or two-dimensional electron gas with low density 
of electrons, Coulomb interaction in the normal region may lead to strong quantum fluctuations which suppress both 
the Andreev conductance and the Josephson proximity effect. Several different kinds of quantum effects are known 
to be relevant in low-dimensional conductors at low temperatures. Quantum corrections to the conductivity of two- 
dimensional dirty conductors grow logarithmically as temperature T decreases and become large at In J*- ~ g (where 
g = (h/e 2 )a is the dimensionless conductance, a is the conductance per square, and t is the elastic scattering time). 
There are two main types of these effects: weak localization corrections [™(|, and interaction-induced corrections 0. 
Other important quantum effects include interaction- induced suppression of the tunneling conductance ("zero- bias 
anomaly" [0J|]), and disorder-induced suppression of superconductive transition temperature T c . These effects 

are of the relative order of g^ 1 In 2 i. e. much stronger than the weak localization and interaction effects. Earlier 
studies of all those effects [p|]7]P,|iO[] employed perturbative diagram technique based on the expansion over semiclassical 
parameter (Epr)^ 1 ~ g ■ This method has an obvious drawback: the number and complexity of diagrams grows 
fast with the order of perturbation theory, which makes its use very difficult in the lowest-temperature region where 
effects are strong. Some combination of the perturbative diagram technique and functional methods was used in 
earlier paper (m), where the first attempt to calculate the effect of long-range Coulomb interaction upon T c was 
made. More advanced functional methods were then developed in the weak localization theory |T^,|l^|, which made 
it possible to average the fcrmionic functional integral over disorder, and reduce it to an "effective" form which 
contains slowly varying diffusive modes only. Those approaches have used either the replica trick [ ^5| or the method 
of supersymmetry |l6|l . Whereas the supersymmetry method was found to be very powerful and convenient for the 
study of single-electron effects, it cannot be used in the cases where quantum corrections due to electron-electron 
interactions are important. The replica method was generalized for the interacting systems by Finkelstein (cf. 
for the review); in particular, he has shown that in dirty films the superconductive T c vanishes at g ~ In 2 (here 
T c o is the bare (BCS) transition temperature). The drawback of the replica method is that it contains an unphysical 
procedure of analytic continuation over the number of replicas n — > 0, and, also, it is difficult to use it for the study 
of non-equilibrium phenomena. 

Long time ago Keldysh flTj] proposed an approach which allows to treat kinetic phenomena in metals with the use, 
both, of the Green function technique, and of the kinetic equation for the distribution function. This approach was 
found to be especially fruitful in the theory of superconductivity, where dynamic equations for the Green functions 
were derived in the dirty limit |f8| (cf. also review articles |f9| , po| ). In the static limit, these equations coincide with 
the Usadel equations [pl| . Keldysh approach often was found to be the most simple and transparent, even for the 
treatment of linear-response problems, since it does not involve tedious analytic continuation procedures, ft is also 
the only known method for the treatment of nonlinear and/or non-equilibrium problems. In some cases nonlinearities 
with respect to both external and fluctuating fields are important, so one needs either to sum up very large number of 
diagrams, or to develop some effective action formalism within the Keldysh approach. Such an approach was recently 
developed, for normal metals, by Kamenev and Andreev p2]] . In many respects we will follow this seminal paper 
in our analysis. A similar approach was also recently developed in p3[ , where Finkelstein's renormalization group 
equations for dirty metals were rederived for the case of short-range electron-electron interaction. For the earlier 
approaches to develop functional integral methods in the Keldysh representation see p4|-p7t . 

In the present paper we will develop the Keldysh functional approach for dirty superconductors, and will use it for 
the study of the Coulomb- interaction effects in the low-temperature Andreev conductance Ga between superconductor 
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and normal metal and in the Josephson proximity coupling Ej via dirty 2D metal. It will be shown that both of the 
above-mentioned effects, interaction-induced suppression of the tunneling DOS, and renormalization of the Cooper- 
channel interaction, contribute considerably into the suppression of Ga and Ej in the low-energy limit. These 
two effects differ in the following sense: specific form of the DOS suppression depends crucially on the long-range 
behaviour of the Coulomb potential (and thus can be varied externally by changing electromagnetic environment), 
whereas renormalization of the Cooper-channel interaction depends on the short-distance Coulomb amplitude only. 
If long-range Coulomb forces are suppressed (i. e. by placing a nearby screening electrode), the DOS suppression 
effect may become weak. In this case the main effect comes from the presence of short-range repulsion in the Cooper 
channel; as a result, both the Andreev conductance Ga{u) and the Josephson proximity coupling Ej(r) exhibit 
anomalous power-law suppression in the infra-red limit, with exponents of the order of <? -1 ^ 2 . In the case of no static 
Coulomb screening, the DOS suppression effect is the strongest one in the asymptotic infra-red limit, it leads to the 
"log-normal" suppression of Ga{u) and Ej{r) as (uj, D/r 2 ) — > 0. The influence of long-range Coulomb interaction on 
the Andreev conductance was treated previously fl28| in a kind of phenomenological circuit theory. In the asymptotic 
limit (l/#)ln 2 (l/W) > 1 our results are in agreement with those of |2S}| and provide a microscopic derivation for 
the effective impedance function used there phenomenologically; in the intermediate region (1/g) ln 2 (l/cjr) ~ 1 an 
additional contribution due to the Cooper-channel repulsion (not treated in pjj ) is shown to be equally important. 

The rest of the paper is organized as follows: in Sec. |l| we describe the formalism of the Kcldysh-type functional 
integral and derive an effective action resulting from disorder averaging; in Sec. Ill we determine the saddle manifold of 
the above action and derive a kind of a nonlinear a- model action formulated in terms of an 8 x 8 matrix Q which depends 
on two times and one spacial coordinate, and of fluctuating order parameter and electromagnetic fields. In section IV 
the basics of perturbation theory and diagram technique for the derived cr-model are presented. In section |y|, in order 
to demonstrate the technique developed, we rederive Finkelstein's renormalization group equations and calculate T c 
suppression for dirty superconductive films. Section VI is devoted to the calculation of the Andreev conductance as a 
function of frequency and/or temperature in the presence of the Coulomb effects; we consider two different geometries 
of N-S contact both in the presence and absence of the DOS suppression effect. In section VII we switch to the 
calculation of the Josephson proximity coupling, for the same two geometrical configurations. Section VIII contains 
discussion and conclusions. Finally, some technical details are presented in two Appendixes. 



II. DERIVATION OF THE EFFECTIVE ACTION IN THE KELDYSH FORM 



A. General procedure 



The Lagrangian of the electron system interacting via electromagnetic field and subject to disorder potential can 
be written as 



where 



and 
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(2.1) 



(2.2) 



(2.3) 



describe electron and field contributions respectively. Here ip a ( r -> t) is a Grassmann spinor field (a = "f , I), <p(r, t) and 
a(r, t) are the electromagnetic potentials, 



(2.4a) 
(2.4b) 



Space integrals in Eq. (2.2) and in the second term of Eq. (2.3) are taken over the system considered while the free 
electromagnetic Lagrangian given by the first term of Eq. ( |2.3| ) contains integration over the whole 3-dimensional 
space, fj, is the chemical potential, and v is the density of states per one spin projection at the Fermi level. A 
short-range interaction in the Cooper channel mediated by electron-electron and clcctron-phonon interactions with 
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momentum transfer of the order of pp is decoupled in the standard way [g9| by the A-ficld. Our superconductive 
coupling constant A coincides with the Finkelstein's definition ]T3] of T c . Generally speaking, one should also introduce 
singlet and triplet coupling constants, T and T2 in Finkelstein's notations, in the diffusion channel. In 2-dimensional 
systems, the correction to them has the relative order of the weak localization effect, (? _1 hi J-, with being the 
relevant energy scale. Therefore, for the study of the Coulomb effects in the Cooper channel and in the tunneling 
density of states, they may be considered constant and can be incorporated into the Fermi-liquid renormalization of 
the parameters of the Lagrangian (2.2). Throughout the paper k = c = 1. 

We are going to construct an effective theory of soft modes in the problem. For this purpose one has to be able to 



take all possible channels into account that is accomplished |16| by introducing a bispinor 



V2 



V 



V'T 

n 



(2.5) 



^ is a vector in the 4-dimensional space 17 which can be considered as the direct product S ® T of the spin tpi) 
and time-reversal (ip, tp*) spaces. The correlations between different time-reversal components of 'J are responsible for 
the quantum correction to the conductivity in the orthogonal case || when the Hamiltonian possesses time-reversal 
symmetry. On the other hand, in studying superconducting phenomena, it is convenient to rearrange components 
of the bispinor in a different manner, separating explicitly the Gor'kov-Nambu p(i| , |3l"[ (ipi,ip*) and spin spaces. 
Finally, one can think of \& as acting in the direct product of the Nambu (N) and time-reversal spaces. These three 
representations are equivalent, ~ S ®T ~ N (g> S ~ N (g)T; we will not specify a certain one and will change between 
them depending on the problem at hand. For later reference, we define the Pauli matrices in the spin, time-r ever sal, 
and Nambu spaces as s,, t% respectively (i — 0,1,2,3). The action of these matrices on the bispinor (2.5) is 
equivalent to multiplication by 4 x 4 matrices; a few examples of them are listed below: 



/o 1 \ 
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(2.6) 



Below we will use the vector VP together with its conjugate vector 5 f+ . These vectors are linearly dependent and 
related by 



1 



where the charge-conjugation matrix 



S = ITy 



In terms of the \& field, the electron Lagrangian C e can be rewritten as 

.„ d (HV - ia> 2 



dr * H 



Of 



2m 



+ U dis (r) + <p + A 



Here 



A = iT y Re A + it x Im A = r + A — t_ A* = \A\[r + e w - T^e~ 10 }, 
is the phase of the order parameter, t± = (t x ± ir y )/2, and the 4x4 matrix S is given by 

2 = t z (g) s = t z (g) s . 



(2.7) 



(2.8) 



(2.9) 



(2.10) 



(2.11) 



Within the Keldysh approach the time-evolution of the system is considered along the Keldysh contour C going 
from t — —00 to t — +00 and then back to —00. At the initial time, t = —00, the system is supposed to be in 
the thermal equilibrium, with the interaction and disorder potential being turned off. The latter are switched on 
adiabatically during the time-evolution. The electromagnetic potentials and a entering the action (2.13) consist of 
fluctuating and external (source) parts: <fi — 4>fi + 0s, a = &fi + a s . The partition function describing the evolution 
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along the contour can be written in terms of the functional integral over fermionic fields ip and ip* and fluctuating 
electromagnetic field potentials <pfi, a/; as 



Z = J D{rMD{<t>f h a } i}D{\}e iS , 



where the action is given by 



S = / Cdt. 



(2.12) 



(2.13) 



External field s (introduced through the source terms in the action) should not be integrated out in the functional 
integral (|2.12j) . If all external fields are classical, i. e. identical for the forward and backward propagation, then the 
evolution along the Keldysh contour brings the system back to the initial state. In this case the partition function is 
automatically normalized to unity, Z = 1. To get a nontrivial result for physical quantities one has to consider the 
generating functional for the source fields having quantum components which are different on the upper and lower 
parts of the contour. We will dis cuss the role of the source terms in section VI B and will operate meanwhile with the 
partition function given by Eq. (2.12). 

Next we divide each dynamical field into two parts residing on the forward and backward branches of the contour 
and labeled by the indices 1 and 2 respectively, and combine them into 2-vectors in the Keldysh space: 



* 2 



ai 
a 2 



A 



Ai 
A 2 



Then the action can be written as S — f_ Ldt with the Lagrangian given by 

L = Ci - £ 2 , Ci = C[ipf,ipi,<f>i,ai,Ai]. 
It is convenient to arrange the components of the Bose-fields (f>, a, A and in the matrix form: 



0i 

02 



a = 



ai 
a 2 



A = 



Ai 
A 2 



0i 
2 



(2.14) 



(2.15) 



(2.16) 



Making use of the Keldysh-space matrices ( 2.16 ) we can rewrite the electron and field parts of the Lagrangian in the 
following concise form: 



L e = / dr* 



, d (SV 



ISl) 



2m 



■/x-Ud is (r) + + A 



(2.17) 



where E and H are expressed in terms of 4> an d a analogously to Eqs. (2.4), and cr, denote the Pauli matrices in the 
Keldysh space. 

Now we are in a position to perform disorder averaging. The Keldysh formalism allows us to average the partition 
function directly utilizing its independence of realization of disorder potential. For the latter we will assume the model 
of a Gaussian ^-correlated white noise with the variance 

([/ dis (r)C/ dis (r')) = ^~ r/) . (2.19) 
Integrating out disorder potential generates a four-fermion term in the action: 

Sdis = ^ / drdtdt>[*+(r,t)aMr,t)] [*+(r, t>,*(r, 0] • ( 2 -20) 



The slow part of the resulting non-local in time action can be decoupled in the standard way (cf. |16|) by the 
Hubbard-Stratonovich matrix field Q: 
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dvdtdt' *+(r, t)Q(r, t, t')cr z *(r, t') - 



8r 



drcftcft' tr Q(r, t, t')Q{v, t' , t) 



(2.21) 



In the Keldysh formalism, Q(r, t, t') is a matrix in the time space as well as an 8 x 8 matrix in the K (g) fi space; it is 
local in the coordinate space once the ^-correlated random potential is considered. Here and in what follows tr(- • •) 
stands for the trace in the K ®Q space whereas the complete operator trace involving integration over space and time 
indices will be denoted by Tr(- • •). 

After the Hubbard-Stratonovich transformation (2.21), the fermionic part of the action becomes quadratic and can 
be written as 



S e = Trf+G" 1 *, 



where the inverse Green function is defined as 

G- 1 



„d_ (sv - csy 

l "~di + 2^ 



The field part of the action takes the form 

E T cr 2 E - H t ct z H 



A 



S f = Tr 
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— Tr A + a z A + — — Tr Q 
A 8t 



(2.22) 



(2.23) 



(2.24) 



As mentioned above, the trace operator Tr(- • •) includes integration over space coordinates; in the first term of 
Eq. (2.24) this integral goes over the whole 3-dimensional space, whereas in the other terms the integral is taken over 
the volume of the system considered. In the present paper we will consider thin metallic films only, so the integral 
in the second and third terms will be effectively 2-dimensional. Gaussian integration over 'J can easily be performed 
resulting in 



S e = — TrlnG" 1 . 
2 



(2.25) 



B. Keldysh rotation 



Among four Keldysh subblocks of the Green function ( 2.23 ) only three appear to be linearly independent [0. To 
simplify its structure, it is convenient to pass to the rotated basis: 



with the unitary matrix L given by pi 



G' = La z GL-\ 



(2.26) 



(2.27) 



After such a rotation G' acquires a triangular form (provided that source terms have no quantum component): 



G' 



G R G K 
G A ' ' 



(2.28) 



where G R( - A ^(t, t') = for t < f (t > f). 

Apart from the Green function, it is also convenient to make a similar transformation for Q: 

Q' = LQL-\ 

and to rotate all 2-vectors defined in Eq. ( [2.14 ) according to 



h - 4>2 



(2.29) 



(2.30) 
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and analogously for a, A and 9 . For the reasons discussed above, (/)[ and <fi' 2 will be referred to as the classical and 
quantum components of the field. The matrices (2.16) will transform according to 



"1 V2 \ — 



(2.31) 



(and analogously for a, A and #), where, following Kamenev and Andreev |22| , we introduced two vertex matrices 

(2.32) 



7 1 = °o, 



7 2 = a x- 



In some cases (e. g., when one considers a uniform superconductor on time scales much longer than the inverse 
gap) it is sufficient to treat the absolute value |A| of the superconductive order parameter as a constant, while 
taking into account fluctuations of its phase. Then, in the rotated basis, the expression for A' can be written as 
A' = |A|[T + e 4ei — T_e~ ie ], or, in terms of the classical, 0[, and quantum, 6' 2 , components: 



A; = \A\[r + e ie ^ - T-e'^] cos 6^ 
A' 2 = \A\[T + e ie 'i +T-e- w 'i}ism6 2 



As a result, the Green function can be written as 



(5V - ia') 2 
% "dt + 2m 



G 



2t 



Q' 



and the action takes the form 



S = --TrlnG'- 1 



Tr 



E'^ovE' — H' T er,,H' 2v-> 
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8r 



(2.33a) 
(2.33b) 



(2.34) 



(2.35) 



The factor 2 difference between the coeffi cients in the terms containing E, H and A in Eqs. ( 2.24| ) and ( |2.35 ) is due to 
the Jacobian of the transformation ( |2.30| ). In what follows we will omit prime at the designation of the Kcldysh-rotated 
fields. This cannot lead to an ambiguity since the original basis will be never used in the subsequent analysis. 



III. (7-MODEL 



In this section we will construct an effective theory that describes low-energy physics of the action (2.35). To start, 
we subject the action to the stationary phase analysis. It is a functional of the matrix field Q and the bosonic fields 
0, a, A, and one has to vary the action with respect to all of them in order to get a set of the saddle point equations. 
First of all we note that quantum components of the bosonic fields are equal to zero in the mean-field approximation. 
Thus, in this section, we will designate their classical components without the subscript "1" for brevity (i. e. (f> = (f>\, 
etc.). Below we will use Q-matrices defined in the energy domain according to the relation 

Q ee , = J J dtdt'e ltt - lt,t 'Qtf. (3.1) 

Varying with respect to Q yields the saddle point equation 

Q(r) = — G(r,r). (3.2) 

In the absence of quantum components, the Green function has a triangular form ( [2.28 ) and so does Q: 



K 



K 



where Q R , Q A and Q K are matrices in the space f2. In the stationary case, the solution Qtv depends on the time 
difference t — t' only, i. e. in the energy domain we have Q e ,e' = 27r5(e — e')Q(e). Varying the action with respect to 
the quantum components 4>2 and a 2 and setting them to zero one obtains the Maxwell equations. In the absence of 
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an external magnetic field and/or voltage drops, the mean-field electromagnetic field vanishes, = a = 0. In order to 
have a closed system of equations, one has to supply Eq. ( |3.2| ) with the gap equation. Varying the action ( 2.35 ) with 
respect to and using Eq. (3.2) we get the selfconsistency equation for the order parameter: 



A = - 



ir\ f de 



(3.4) 



where the 4-dimensional space fl is defined in Sec. II A 



To clarify the structure of the saddle point given by Eqs. ( |3.2|) and (3.4), consider first the case of nonsuperconducting 
(A = 0) metal. Then it is easy to check that a diagonal in the energy space matrix 



A(e) = A (e)3, 



where 



A (e) = 



1 2F(e) 
-1 



(3.5) 



(3.6) 



K 



is a solution of Eq. ( |3.2| ). The 4x4 matrix function F(e) introduced in Eq. (3.6) has the meaning of a generalized 
distribution function. In the steady state it reduces to the single scalar function F(e) = 1 — 2/(e) where /(e) is the 
usual Fermi distribution function. The form of /(e) is not fixed by the saddle point equation since A(e) with arbitrary 
F(e) satisfies Eq. ([T^). The reason for this ambiguity is that any distribution function is allowed in the absence of 
interaction. To bring the system into the equilibrium with 



F oq (e) = tanh 



2T 



(3.7) 



either electron-electron or electron-phonon inelastic interactions must be taken into account. The Keldysh formalism 
is suitable for the study of noncquilibrium problems as well. In this case there is an externally controlled differen ce o f 
temperature and/or chemical potential across the system, and the function F(e, r) should be obtained from Eq. (3.2) 
with the proper bou ndary conditions. 

The solution (|3.5| ) captures the eigenvalue structure of a generic saddle point. All fluctuations of Q that alter the 
eigenvalues ±1 are massive. The massless modes share the eigenvalue structure of A and can be obtained from it by 
the following transformation: 



= U- X AU, 



(3.8) 



where U is some rotation matrix which acts in the 8x8 space K (g> fl as well as in the time domain. According to 
Eq. (3.8), the field Q satisfies the nonlinear constraint 



Q 2 = 1 (3.9) 

at the saddle point manifold (SPM). Together with Eq. (|3.3|), this suggests the following parametrization of the 
Keldysh block: 



Q F - FQ , 



(3.10) 



where, again, F has the meaning of a generalized distribution function. 

Now let us turn to the case of a uniform bulk superconductivity. Here it is convenient to chose a representation 
of the space VI as a direct product of the Nambu and spin spaces, f2 = N ® S; in this notations S = t z . The 



superconducting saddle point solution, Qs, has the form (3.3) with fll 



±- 



1 



V(e±iO) 2 -|A| 2 



e A 

-A* -e 



(3.11) 



N 



Taken at the saddle point manifold, the matrix Q is equivalent to the Larkin-Ovchinnikoy [lq| quasiclassical Green 
function g. The mean-field value of the order parameter can be obtained form Eq. (3.4). Substituting Q K from 
Eqs. (|3T0| ) and ( |3~TT| ) we obtain the standard BCS gap equation (negative A corresponds to attraction between 
electrons) 
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A = -AA 



de e 
— — tanh — — . 

^e 2 - |A| 2 2T 



(3.12) 



The superconducting saddle point, Qs, belongs to the metallic SPM given by Eq. (3.8). However, in the presence 



of superconductivity, some excitations on the metallic SPM, having been massless, acquire a gap proportional to A. 
A detailed discussion of the hierarchy of gaps in the cr-model for N-S systems can be found in the review 32 1. 

The next step in the derivation of the er-model is to consider fluctuations of the SPM and to perform the gradient 
expansion of the action ( [2.35 ). Such a procedure is justified in the dirty limit, At -C 1, which will be implied 
from now on. It is equivalent to the replacement of the full Eilenberger j3j| equations in the conventional theory of 
superconductivity to their approximation proposed by Usadel [£1[ . 

As was recently suggested in Ref . |22| , in studying electric field fluctuations it is convenient to single out the gauge 
degrees of freedom in Q by the transformation 



Qw 



= e mm ( 



we 



-iK(t')3 



(3.13) 



where K = Kij 1 is related to the doublet K = (K\, K2) T , in analogy with the field cf>. After the transformation 



(3.13), the action can still be written in the form (2.35), with the Green function being substituted by 



G~ 



(SV - iAf 

2m 



2r 



Q + $ + A K , 



with 



A = a - VK, 

5 = % - d t K, 

Expanding TrlnG -1 in the standard way MJ22J, we obtain the following effective action 



S = S a + 2^Tr$ T cr x $ + Tr 



E T cr T E - H T cr r H 



Aire 2 



yTrA+^A, 



* + A x )g 



where S a is the cr-model action for the matrix field Q, 

S a = ^Tr\D(dQ) 2 + 4i(iZd t 

o 

Here D is the diffusion coefficient and d denotes a long covariant derivative, 

3X ee VA - i [SA, X] . 



(3.14) 

(3.15a) 
(3.15b) 
(3.15c) 

(3.16) 
(3.17) 
(3.18) 



Derivation of the effective action (3.16) that describes interacting disordered normal/superconducting electron liquid 
is the main result of this section. The model is formulated in terms of the interacting matter field Q subject to the 
nonlinear constraint Q 2 — 1, electromagnetic fields (j>, a, and the pairing potential A. At the present stage, the 



phase K introduced in Eq. (3.13) is left unspecified. It will be expressed in terms of the electromagnetic potentials in 
section [V C . 

x 8 matrix in the K H space. In what follows we will assume that all interactions are 

)S 



In Eq. ( |3.17| ), Q is an 

spin-independent. Then Q is proportional to the unit matrix in the spin space, and the 4-dimensional space £1 = 
collapses into the 2-dimensional Nambu space. As a result, the theory will be formulate d in t erms of the 4x4 matrices 
Qtt> acting in the K ® N space. The corresponding action can be obtained from Eq. ( 3.17 ) by taking the trace over 
the redundant spin space: 



S a = — Tr 



D{dQ) 2 + 4i(ir z d t + $ + A K )Q 



where the operator d is given by Eq. ( 3.1 8| ) with S = t z . 

Varying the action Eq. ( 3.19| ) with respect to Q with the constraint Q 2 — 1 yields the equation 



Dd(QdQ) + i [ir z d t + $ + A K , Q] = 0, 



(3.19) 



(3.20) 



which (for K — 0) coincides with the dynamical Usadel equation |l8| ]. In the absence of superconductive coupling, 
A = 0, Q is proportional to r 2 : Q = Qkat z , and our action (3. IS) reduces to the Kamenev-Andreev action p2j for 
the field Qka- 



9 



IV. PERTURBATION THEORY 



A. Free metallic diffusons and Cooperons 



In this section we will show how a systematic perturbative expansion of the er-model ( 3.1 6| ) can be devel oped . 
Keeping in mind further application to N-S devices with relatively weak proximity- induced coupling (sections VI, VII), 
we will consider fluctuations near the metallic saddle point (3.5). For this purpose, it is convenient to parametrize 
the rotation matrix U in Eq. (3.£) in terms of another matrix W subject to the linear constraint 

WA + AW = 0. (4.1) 

Such a parametrization is not unique and a number of them are widely used in literature (see, e. g., [fl6|| ). 

In the Keldysh formalism, the saddle point ([si]) is not diagonal and, consequently, the solution of Eq. (hi) for 
W explicitly depends on the distribution function F(e). As a result, even for a noninteracting system, intermediate 
expressions for the Cooperon and diffuson propagators would depend on the particle distribution. To surmount such 
an unphysical complication, we note that at the saddle point ( p\3| ), Q can be diagonalized in the Keldysh space by a 
nonunitary 2x2 matrix 



1 F 

-1 



(4.2) 



K 



that separates the distribution function F from the retarded and advanced blocks which are determined by the spectral 
properties only: 



Q = u 



Q R 
Q A 



(4.3) 



The function F in Eq. (h2) is the stationary fcrmionic distribution function and we will assume that the system is in 
thermal equilibrium, s o th at F(e) is given by Eq. (3.7). 

The decomposition (O) suggest to pass from the initial Q-representation to a new variable, Q, defined as 



In terms of the new variable Q, the a- model action ( |3.16 ) acquires the form 



S a = — Tr 



d <-> «-» 
D{dQf +Ai(iT z — +u<Pu + uA K u)Q 



with the modified definition of the long derivative: 

dX = VX -i[r z uAu,X], 



(4.4) 



(4.5) 



(4.6) 



Note that the matrix u couples to the interaction terms only. For the noninteracting case, the distribution function 
drops from the Q-action. 

At the metallic saddle point for the action (4.5), $ = A = A = and Q is diagonal, Q = uK^t z u = <j z t z . Gapless 
fluctuations of Q can then be parametrized as 



Q=U- 1 a z r z U, 

where U is a unitary matrix, for which we adopt the exponential parametrization, 

U = e w / 2 , 

in terms of the matrix W which anticommutes with <j z t z , 

{W,a z r z } = 0. 
An explicit expression for Q in terms of W reads 

Q = e- w / 2 a zTz e w ' 2 = a z r z (l + W + W 2 /2 + ...). 



(4.7) 
(4.8) 
(4.9) 
(4.10) 
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The linear matrix constraint (|4Jj) can be resolved by introducing eight scalar variables, Wi and Wi with i = 0, x, y, z, 

as 

\ Hl„ _L 111 T 111 1- —1— 1/1 T" / \ / 



W + U> Z T Z WjT,. + WyTy 



K 



The diagonal (offdiagonal) in the Nambu space excitations, Wi and Wi with i = 0, z (i = x,y), correspond to diffusion 
(Cooper) modes. Extracting quadratic in W part from the cr-model action ( [4.5D , 

i^ 2 ' [W] = ^ Tr [(Dg 2 - 2i<7,E) W(q)W(-q)] , (4.12) 

we obtain the following correlators for diffusons: 

(w i (q;ei,e 2 )w i (-q;e3,e 4 )) = v "' v ^^ 1 77T~"^~T \ ~ ' « = 0, z, (4.13) 

and Cooperons: 

(mi(q;ei,e 2 )wi(-q; £3,64)) 

i = x,y. (4-14) 

I |7t|-,»I( , — /• , Ml,., — I 

(Wi(q; ei, e 2 )uJi(-q; £3, e 4 )) 



1 


(27T) 2 5( £l 


- £4)^(^2 - 


e 3 ) 


nv 


Dq 2 - 


- i(ei - e 2 ) 




1 


(27r) 2 5( £l 


- £4)<^(£2 - 




7T// 


Dg 2 - 


- *(£i + £2) 




1 


(27r) 2 5( £l 


- £ 4 )^(£ 2 - 


e 3 ) 



7T^ Dq 2 + j(£l + £ 2 ) 



B. Diagrammatic technique 



Expression ( 4.1 0| ) provides a regular way for perturbative expansion near the metallic saddle point, Q = <j z t z . Its 



basic elements are given by the free correlators ( |4.13|) and ( 1.14 ) corresponding to soft diffusion and Cooper modes 



Since W ee ' is a matrix in the energy space, its correlators are represented diagrammatically by two parallel lines, 
each of them carrying one energy index, see Fig. 0a. Expanding the action ( |l.5| ) over W generates nonlinear vertices 
describing interaction between diffusion and Cooper modes. The resulting diagrammatics looks very similar to the 
standard cross technique for dirty metals fjij where soft modes are constructed from two Green functions averaged 
over disorder. Note however that in the present case both diffusons and Cooperons are depicted in the same manner, 
with arrows pointing in the opposite directions on two lines of the propagator. Such a convention, though being 



unusual for Cooperons, is consistent with the definition of the Fourier-transformed variables (3.1). Diffusons and 
Cooperons can be distinguished by their structure in the Nambu (and Keldysh) space. 

£ 2 



(a) W £lEl , 
£1 

(b) A -<^ Q (c) Q X Q 

FIG. 1. Basic elements of the diagrammatic technique: a) diffusion/Cooper propagator; b) Q - A interaction vertex; c) 
effective interaction S\[Q] after eliminating A field. 

The normalization of the functional integral, Z = 1, manifests itself in cancellation of closed loops in the perturbation 
theory. In the Keldysh formalism, such a cancellation is related to the integral over the internal energy (let it be E) 
of a closed loop. Indeed, all propagators along the loop have poles in one (upper or lower) half-plane of the complex 
variable E. Therefore, integrating over E yields to cancellation of the corresponding diagram. 

For future references we present here the contraction rule for averaging over W: 

1 f dqde 1 de 2 J tr(AB - Acr z Ba z + At z Bt z - Aa z T z Ba z r z )(Dq 2 + i{e x - e 2 )a z ) 



, tx(AB + Ao z Bo z - At z Bt z - Aa z r z Ba z r z ){Dq 2 + i(ei + e 2 )a z ) | ; [ { _ . 



(ZV) 2 + (£i + £ 2 ) 2 
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where A — A(q,, £1,62) and B = B(— q, £2,61). The first (second) term corresponds to diffusons (Cooperons). For 
Dq 2 3> £i, £2, this expression can be simplified as 

(Tr AW ■ Tr BW) = - J- / ^ ~ A °f* B °^ , (4.I6) 

ttu J (2Tr) d+z Dq 1 



Apart from the matter field Q (or, equivalently, W), the cr-model action ( p.lq ) with So- given by Eq. ( |3 . 1 9[) contains 
the electromagnetic potentials and the field A. The former will be considered in the next subsection while the latter 
can be easily eliminated by Gaussian integration. Note that the resulting expression does not depend on the Coulomb 
phase K entering Ak since it can be "gauged away" by the shift of the phase of the integration variable A. In other 



words, Ak in Eq. fl3.19|) can be substituted by A provided that the order parameter field is to be integrated out. The 



resulting contribution to the action reads 

7T 2 i/A 



drdt tenQtt ■ tona x Q tt = / drdttra x [Q 2 - (r 2 Q) 2 ], (4.17) 



and is shown diagrammatically in Fig. [l|c. Note that this term is conveniently expressed in the original Q- 
representation. 

The perturbative expansion near the metallic saddle point is justified in the case of repulsive (A > 0) interaction in 
the Cooper channel or when the proximity-induced superconducting coherence is weak enough. Otherwise, deviation 
from the metallic point is large and one should use the solution to the Usadel equation as a starting point for the 
perturbative analysis (cf. a detailed discussion in Ref. f32j). 



C. Electromagnetic field fluctuations 



In the previous section we have sketched the basic rules of the diagrammatic perturbation theory in the absence 
of the Coulomb interaction. Here we will discuss how fluctuations of the electromagnetic field can be incorporated 
into the formalism. As was shown recently by Kamenev and Andreev p2|, a certain choice of the Coulomb phase 



K introduced in Eq. (3.13) results in a significant simplification of the theory. Such a choice essentially depends 
on the position of the noninteracting saddle point on the SPM and is quite different for the metallic (|]^) and 
superconductive (3.11) saddle points. For the metallic (nonsuperconductive) case, we choose, following Ref. p2[ , K 
to be a linear functional of cj) and a and require the vanishing of the term linear both in W and A in the Usadel 



equation ( [3.2C ), or, equivalently, in the cr-model action (3.19). The resulting equation reads 



£(VA-A VAAo) + [A ,$] = 0, 



(4.18) 



where the matrix Ao is introduced in Eq. (|3.6| ) . Eq. ( 4.18 ) is to be used to express K in terms of the electromagnetic 
field potentials <j) and a, the corresponding relation having the form (in this section we share most of notations of 
Ref. @) 



V 1 (uj,q)K{q 1 uj) = IT W V(q,^) + iDa x qa(q,u>). 



Here 







Dq 2 + iuj 



IT 



Dq 2 - iuj ~2iujB u 

1 1 -1 

1 25 w 



and 



B,, 



coth - 



2T 



(4.19) 

(4.20) 
(4.21) 

(4.22) 



is the equilibrium bosonic distribution function. 

So far our analysis holds for any gauge and any geometry of the sample. From now on we chose the gauge a = 
neglecting relativistic effects due to the magnetic field fluctuations. Also we restrict ourselves to the consideration of 
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2-dimensional systems. Then one has to integrate out-of-plane degrees of freedom in the electromagnetic field action 
S em = TrE T cr 2: E/(47re 2 ). The result depends on the presence or absence of conducting electrodes that screen the 
long-range Coulomb interaction (for the sake of simplicity, we set the dielectric permeability of the medium to unity) . 
For a single plane, one obtains 

dq 



dt 



(2tt) 



•tr^-qKV^tq), 



where Vo(q) is the 2D Coulomb interaction potential, 

t r / \ f dk z Aire 2 



2ne 2 



(4.23) 



(4.24) 



If, for example, there is a metallic gate at a distance b from the 2D plane, then Vo(q) is screened in the long- wavelength 
limit and we have instead 

2ttp 2 

V Q s "( q ) = _ (1 - e - 2b «) . (4.25) 



Collecting the terms in the action bilinear in and K and making use of the relation ( 4.19 ) we get an effective 
action for the electromagnetic field propagation in the disordered metal : 

S%tt = Tr 1 T V~^, (4.26) 

where V has the meaning of a dynamically screened Coulomb interaction in the RPA approximation, 

V(q,L0) = (a x V Q -\q) + Poiq^))- 1 , (4.27) 

where Po(q,u) is the bare density-density correlator. The matrix V(q,uj) has the structure of a bosonic propagator 
in the Keldysh space |p2| : 

V K {q,u) V R (q,u) 
V A (q,uj) 



V(q,u) = 



(4.28) 



with 



2vDq 2 



V R > A (q,u;)=(v -\q)+ , , 
V K (q,w) = B u (V R (q,Lj) ~ V A {q,u)) . 



An additional factor 2 in Eq. (4.29a) compared to that in the Kamenev-Andreev paper 
they considered spinless electrons. 

Eq. (4.26) determines the bare propagator of the electromagnetic field: 

(<fo(q,w)</>j(-q,-w)) = -Vy(g,w). 



(4.29a) 
(4.29b) 

is related to the fact that 



(4.30) 



The propagator of the field K can be obtained from Eq. (4.3C) with the help of the relation (4.19) and has the form 



(Kifa^Kji-q, -w)) = -Vij(q,u)), 



where the matrix V has the same structure as V, Eq. (4.28), with 



V 



R,A 



few) = - 



1 



{Dq 2 =F iuj) 



V -\q) + 



2vDq 2 
Dq 2 =p ito 



(4.31) 

(4.32a) 
(4.32b) 



V K (q,w)=B u) (V R {q,u)-V A (q,u>)). 

We should note, that in Eqs. ( [4.29a ), ( [4.32a] ) possible influence of superconductive pairing upon the dynamic screening 
of the Coulomb interaction is neglected; this is safe since we will consider N-S systems in the limit of weak tunneling 
only, so the 2D metal is slightly perturbed by superconductivity. 

In a superconductor, the choice of an optimal Coulomb phase K valid in the whole energy range is a complicated 
task. However, it had been shown in Ref. |35|] that in the deep subgap limit (e ^ A) the effect of the electric potential 
on the quasiclassical Green function Q is small in the parameter e/A and hence K = 0. This result will be used 
below. 
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V. RENORMALIZATION OF THE INTERACTION IN THE COOPER CHANNEL 



A. Renormalization Group procedure 

In this section we will show how to construct a procedure of successive eliminating of high-frequency and high- 
momentum fluctuations of all interaction modes in the problem: the matrix field Q(q, e, e'), the order parameter 
field A(q, uj), and the electric potential </>(q, uj). Elimination of high-energy modes in a dirty 2D metal results in 
logarithmic corrections to the parameters entering the action and governing dynamics of the retained slow modes. 
This procedure is known as the Renormalization Group (RG) method. We will closely follow Finkelstein's approach 
to RG construction jy], with the necessary modifications due to the presence of the Keldysh space instead of replicas. 

Each elementary RG step consists in elimination of degrees of freedom in the energy shell from fi* to Q, where 
O is the current value of the running ultra-violet cutoff in the problem. Correspondingly, all fluctuating fields are 
decomposed into fast (denoted by a prime) and slow (denoted by a tilde) parts. For the fields A, </> and K such a 

representation is trivial: A = A + A', etc., while for the field Q it must be consistent with the constraint Q 2 = 1. 



To achieve this, we decompose the rotation matrix U in Eq. (4/7) into the product of a fast, U' = exp(W'/2), and a 
slow, U, part, so that 



Q=U- 1 Q , U, 



(5.1) 



where the fast Q' is expressed in terms of W' according to Eq. ( 4.1 0| ) . The slow matrix U differs from the unit matrix 
only if all its arguments are smaller than the new cutoff f2*: 



U(q,e,c') = l, if Dq 2 or |e| or \e'\ > 



(5.2) 



On the other hand, the fast W is nonzero if at least one of its arguments belongs to the energy shell (£7*, f2): 

W(q, e, e') ^ 0, if fi* < Dq 2 or |e| or |e'| < a (5.3) 

After integration over fast variables, f2* becomes a new cutoff and the whole procedure should be successively repeated. 

In a 2D dirty metal, integrating out fast degrees of freedom results in a relative correction oc ln(Q/Q*) = — £ to 
the parameters of the effective action (3.16), where £ is a logarithmic variable defined as 



C = ln 



(5.4) 



The RG near the metallic saddle point a z T z is justified provided that the Cooper-channel coupling constant A <C 1 
while the dimensionless conductance of the metal </> 1. The latter is defined as 



g = 2vD = 



(5.5) 



with a = -Rq 1 being the conductance per square (in conventional units, a — (e 2 /K)g). Logarithmic corrections to the 
conductance become large at the localization scale C ~ 9- The same is true for the coupling constants T and T2 omitted 
in the derivation of the action ( 3.16 ), cf. discussion in Sec. [I A . On the contrary, corrections to A become of the 
relative order of unity at much shorter scale, at £ ~ ^fg. Therefore it is possible, at large g, to neglect renormalization 
of the conductance and consider g as a constant. 



B. BCS correction to A 



First of all we show how to obtain the standard BCS renormalization of the Cooper-channel interaction constant 
A in the present formalism. The correction originates from the term Saq = — 7wTr AQ in the a- model action ( 3.19| ) 

after eliminating high-frequency fluctuations of the field Q (here we substitute by A as was explained in section 
IVB). Passing to the rotated Q-representation (4.4), expanding to the first power in W" according to Eq. (4.10) and 
setting U to unity one obtains the relevant interaction vertex 



-iTv Tr uAua z T z W' . 



(5.6) 
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After averaging over the fast W', this term will generate the following correction to the action: 



A5 = -(S?, 
8+C0/2 



(5.7) 



>- -®,-q 



8-C0/2 

FIG. 2. BCS correction to A. 

The corresponding diagram is shown in Fig. pi The only fast variable is the internal energy of the (WW") 

propagator. Employing then the contraction rule ( j4.15| ) and the relation {A,t z } = 0, we conclude that the diffuson 
contribution vanishes identically while the Cooperon pairing yields (to logarithmic accuracy) 



AS 



■kv 

T 



dq dw de 1 

(2tt) 4 7 



trA(q,w)A(-q, -w)A (e) 



(5.8) 



where the e- integration is perfor med over the energy shell fi* < |e| < Q, while q and uj are restricted to the domain 
Dq 2 , \u>\ < 57*. On deriving Eq. (5^8) we have used the relation Aq = ua z u. Now using the definition of the matrices 



A and Ao and omitting the tilde sign over the designation of the slow component of A, we transform the trace in the 
above equation to the form 



tr(r+A 4 - r_A*)(r+A j - r_A*) 7 ¥A (e) = -AF{e){A\A 2 + A,A*). 



(5.9) 



As a result, Eq. (5.8) can be represented as 



— » — » f d 

AS = 2vTiA+a x A- — tanh— . (5.10) 
The temperature T thus determines the infrared cutoff of the RG procedure. At larger scales, 51 3> T, one has 

AS = 2u In — Tr A + <r x A. (5.11) 



Comparing with the last term in the action (ft.lq 1 ), we conclude that the correction (5.11) renormalizes 1/A: 



1. 



(5.12) 



For a superfluid Fermi-liquid, the RG equation (5.12) was derived in j pq ]; in that case only the Cooper channel was 
relevant, therefore the RG approach was equivalent to a simple summation of the standard BCS-theory ladder. 

Coulomb interaction in dirty superconductors contributes both to the Cooper and to the density-density channels. 
We consider here the range of parameters where ln(l/Qr) <C g that makes it possible to neglect the effect of the 
Cooper channel upon the conductance g; however the effect of the density-density channel upon the Cooper one has 
the relative order of g^ 1 In 2 (1 /fir) and thus should be taken into account. This effect can be described in the form 
of an integral equation for the energy-dependent Cooper attraction A(-E), as was done by Aleiner and Altshuler |37f | 
(the same kind of equation was derived in for another but similar problem). For our purpose it will be more 
convenient to treat the same effect within the RG procedure, as described in the next subsection. 



C. Coulomb correction to the Cooper-channel interaction 

In this subsection we calculate the Coulomb-induced correction to the coupling constant A. It appears as a result 



of eliminating high-momentum fluctuations of the electric field. According to Eq. (4.5), the electric field couples to 
the matter field Q by the following terms: 



Sint = -ttvTt 



u{4> ~ d t K)uQ + ^Dt z uWKu[Q,WQ] + 1 -D[t z uVKu, Q]f 



— °int 



Ob 
Mr 



(5.13) 
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Here we utilize the relation ( 4.1 9| ) connecting the phase K to the field (j>. The interaction vertices Sf ut ana ^ ^fnt are 
linear in 0, while the vertex Sf nt is quadratic. Then the result of averaging over fast variables can be written as 



As = -([sr nt + sfj 2 > + (sf nt >. 



(5.14) 






(C) 

FIG. 3. Diagrams for the Coulomb correction to A: a) ((S^) 2 ), wavy line denotes the correlator {4>i(<\, ui)(f)j(—q, — lo)); b) 
{(Sf nt ) 2 ), zigzag line stands for the correlator (-K;(q, u))K 3 ■(— q, — uj)); c) (Si C nt ). 



Consider first the average ((S 1 ? 



The relevant (i. e. logarithmic) contribution can be written in the form 



((Sf nt ) 2 ) = 7rV( [Tr Uu$ - d t K')uU~ l a z T z W] 



(5.15) 

where the pairings are shown in Fig. |3^,. According to this diagram, all energies are coupled to slow variables, and 
there is only one fast variable, the internal momentum q running over (4>'<j>') and (W'W). Therefore, to logarithmic 
accuracy we may consider Dq 2 3> u>. In this limit, one can neglect the term iojK compared to cf> as follows from 
Eq. ( 4.19 ) and use the universal large-*? asymptotics of the screened Coulomb interaction ( 4.28| ), 



»(q J w)(/) i (-q,-w)) 



4i/ 



(5.16) 



which is independent of the details of the bare potential Vo (q) . The averaging over W' can performed with the help 
of Eq. ( 4.16|) , making use of the inequality Dq 2 » e%, e 3 . As a result, one obtains 



((SLY) = 



dr dei de% de^ de± du> dq 1 



Dq' 



^ tr(AiBj - Aia z T z Bja z T z ), 



(5.17) 



where the q-integration is taken over the fast energy shell f2* < Dq 2 < il, the matrices A{ and Bj are given by 



B, 



<=3 



(5.18a) 
(5.18b) 



and all slow matrices U are taken at a coincident spacial point r. Performing integration over e\ and £3 in the first 
term under the trace in Eq. (5.17), we obtain 



—^-Y j txA i B j = 2txa x {uQu) i 



C4 + LJ,€2 ) 



(5.19) 



where a slow Q is defined as Q = Li a z T z Li. In the second term under the trace in Eq. (5.17), the matrices o z t z 
cancel with those in Eq. ( j.lq 1 ) , integration over e\ and £3 is equivalent to multiplication of U~ 1 and U that gives the 
unit matrix, and tracing with a x yields zero. 
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The resulting expression takes a simple form in the initial Q-representation in the time domain. Integrating over 
q, one obtains (omitting the tilde sign) 



iv . 



((S?ntr> = -77 In 



4g Si* 



drdttr a x Q 2 tt iv). 



(5.20) 



The correction to the action (i/2)((5'int) ) ^ s the oru y one m the standard gauge with K = 0. Physical quantities 
should not depend on the choice of K, which is a kind of a gauge transformation upon the matrix Q (cf. Eq. (3.13)). 
For this reason, it would be enough to use Eq. ( 5.20 ) to find the renormalization of the coupling constants in the 
effective action. Below we present, however, the calculation of other terms in Eq. ( 5.14 ) related to K, in order to 
show that they really cancel each other. 
The diagram for 



((Sf nt ) 2 ) = 7rV£> 2 ( [Tr Ut z uVK'uU- 1 VW'Y 



(5.21) 



is shown in Fig. ||b. It looks similar to the one in Fig. ^a, with the wavy line being replaced by the zigzag line denoting 
(K'K 1 ) correlator. Repeating the steps that lead to Eq. ( pM7| ) and using the large-g asymptotics of ( 4.32 ) we obtain 

(5.22) 
with 



, lQ b n2\ i7T f dvde l de 2 de 3 de A dudq_ 1 ^ 

((Sint) ) = "J / T^w j^Wi ~ Mi<J z T z Nj(J z T z ) 



Mi = U eiea T z u ea j l u ea - u U e i 



£2 — ",£3 ' 



Integrating over q, t\ and £3 as described above, we get 



((SLY) 



iv j f2 
4# Si* 



dr^trcr x (r z (9 tt (r)) 2 



(5.23a) 
(5.23b) 

(5.24) 



In the same manner it can be shown that the average (S^t'^int) = vanishes. 

In calculating the average (S^ nt ) shown diagrammatically in Fig. |]c, all Q's may be considered slow. The analytical 
expression reads 



(5P t ) = -^(Tr(r,QV^ 



'\2 



(5.25) 



where the fast momentum runs over the (K'K 1 } propagator. Calculating the corresponding logarithmic integral, one 
obtains 



v £1 f 

(Sint) = ~ ln ft~ J dlC dt tr a * ( T zQt 



«) 5 



(5.26) 



Substituting Eqs. (|5.20|) , (|5.24|) and ( |5.26|) into Eq. (|5.14|) , we see that the contributions from the vertices describing 
interaction with the field K cancel each other, and the resulting expression is given by 

v si f v si f 

A5 = ^ln— drdttra x Ql(r) = —ln— / dr dt V tr a x TjQ tt (r) ■ t 3 Qu(t). (5.27) 
85 J 16g si* J ^ 



Comparing with Eq. (4.17), we see that the terms with j = x,y renormalize the Cooper channel coupling A, whereas 
the terms with j = 0,z contribute to the couplings in the diffusion channel. The latter, T and T2, are not taken 
into account since corrections to them are of the relative order of g~ x ln(l/ilr). As a result, we get the Coulomb 
contribution to the RG equation for A: 



d( 4ttV 



(5.28) 



which coincides with the Finkelstein's result 1 1 1 13] in the limit T, T2 —>■ 0. 
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D. Solution of the RG equation for A and shift of T c 



Combining Eqs. (5.12) and (5.28), we arrive at the complete renormalization group equation for the Cooper-channel 
interaction constant: 



— - -A 2 — !— 



(5.29) 



In the high-energy range 1/r < E < Ep the second term of Eq. ( |5.29] ) is absent. We will use the solution of Eq. (5.12) 
at E sa r" 1 as an initial condition for t he ful l equation ( 5.29 ). 

The renormalization group equation ( |5.29 ) possesses two fixed points, ±A S , where 



(5.30) 



The stable fixed point, +X g , is the limiting point of the RG flow in the metallic region. A trajectory reaches its 
asymptotic value A g at the scale £ ~ Jg. The unstable fixed point, — A 9 , separates the regions of metallic (A > — X g ) 
and superconducting (A < — X g ) states. The solution of Eq. ( 5.29 ) is given by 



A(C) 



Ao + A s tanh X g ( 

1 + — tanh X g ( 
A s 



(5.31) 



where Ao is the bare value of the interaction constant in the Cooper channel defined at the energy scale r . 

To study the superconductor-metal transition we will consider here the case of an attractive interaction, Ao < 0. 
The superconducting transition temperature T c is determined by the position of the pole in A(£), with ( = In j, — . In 
the clean system (g — > oo), 

(5.32) 



T c0 r = exp(- — 



For a finite X g , the critical value of £ can be easily obtained from Eq. (5.31) and is given by 



Cc 1 JX \ + X g 



2A„ 



|A |-A S 



Consequently, we get for T c : 



T r r = 



|A |-A 3 
|Ao|+A 9 



(5.33) 



(5.34) 



Substitution of A in terms of T c q with the help of Eq. (5.32) leads to the final result for T c suppression by disorder, 
which coincides (within our accuracy) with that of pi: 



T c t = 



1 + 



T c0 t 
1 ln^- 

27T^g T c0 r , 



TVs 



(5.35) 



Evaluating expression (5.35) for g> 1, we obtain a perturbative reduction of the transition temperature |9j,|l 



In- 



T, 



ln J 



rO 



12^ 2 «? T c0 T' 



(5.36) 



valid at large conductances. The critical temperature becomes zero and the superconducting transition disappears at 
the critical value of the dimcnsionlcss conductance 



g c = 



2tt ln T c0 t 



(5.37) 



Note once again, that the result (5.35) for T c (g) dependence is obtained neglecting weak- localization corrections 
to the conductance g, as well as thermal, quantum and mesoscopic fluctuations. This is correct provided that the 
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renormalized conductance at the scale T c is still greater than unity: g* = g — {a/-K 2 ) In 3> 1, where the constant 
a is given by the sum of the usual weak localization and interaction corrections, and is equal either to 1 or 1/4, 
depending on the absence or presence of the spin-orbit interaction. 

In the above derivation we neglected spin-dependent interactions (T^ in the notations of ]T^|) which may change 
the numerical coefficient in Eq. (5.28). On the other hand, it was explained in p"2| that strong spin-orbit scattering 
eliminates possible effect of T2 upon T c (g) dependence. 



VI. ANDREEV CONDUCTANCE 

A. Tunneling term in the action 

In the previous sections, we considered a uniform 2D system. In principle, spacial inhomogcncitics in the local 
system's characteristics such as the conductance and/or the Cooper channel interaction A can be easily incorporated 
into the er-model action ( |3.16 ) by a spacial dependence of the parameters of the cr-model. Thus, the action (3.16) is 



suitable for a description of N-S interfaces or interfaces between metals with different conductances. However since 



the solutions of the Usadel equation (3.20) are continuous, only interfaces with perfect transmission T — 1 can be 
described in such a manner. In order to be able to deal with the interfaces of arbitrary transparencies, one has to 
introduce a boundary term into the action. Below in this paper we consider the case of low-transparent interfaces 
which can be described by means of the tunneling Hamiltonian approximation. Then the boundary term in the action 
can be derived in the second order over the tunneling Hamiltonian and reads 

S, = f 7 Tr r g«Q( 2 ). (6.1) 

Here and refer to different sides of the interface boundary T; the notation Trr means that the space integral 
is taken over the interface surface, and 7 is the (dimensionless) normal-state tunneling conductance per unit area of 
the boundary. In the expression ( |3.l| ) the trace over the spin space has already been performed. 

Variation of the total action 5 to t = S& + Sy, with So- given by Eq. ( |3.19 ), with respect to Q, reproduces the matrix 



Usadel equations for Q*- 1 -* and together with the corresponding boundary conditions |3£]j3^] (cf. similar derivation 
in®): 

gi QWV^=g 2 Q^V^ = |[g< 2 \Q«], (6.2) 

where V_l stands for the gradient along the normal to the interface directed from the medium (1) to the medium (2) 
(Eq. (pj) is written in the form assuming the absence of magnetic field) . It amounts to a straightforward calculation 



to show that the action (6.1) leads to the following expression for the bare normal-state tunneling conductance ax- 



a T = —Ay, (6.3) 
h 

where A is the area of the tunnel junction. We omit here such a calculation since it is fairly similar to the calculation 
of the Andreev subgap conductance presented in the next subsection. Similarly, we will not dwell upon Coulomb 
interaction-induced corrections to the tunneling conductanc e ]7| , [4l| ] which can be derived from the action ( |6.1[ ) by 
taking into account fluctuations of the fields Q, <j) and K (cf. ]22|). An analogous calculation of the interaction effects 
in the Andreev conductance is one of our main subjects below. 



B. Andreev conductance in the effective action formalism 



In this subsection we redcrive, within the effective action formalism, the well-known results for the subgap Andreev 
conductance Ga between a superconductor and a dirty normal metal (cf. e. g. |Q). We start from the simplest 
situation when Ga does not depend on voltage and/or frequency. First of all we show that if the effective action 
contains the following term 

Sa = ^G a Tt(Q s A) 2 , (6.4) 
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then G a is indeed the dimensionless subgap conductance. In the second step, we prove that the term of the form 
( |6.4[ ) is generated in the second order of expansion over S 1 . 

To begin with, we note that at low energies, e <C A, the superconductive matrix Green function Qs does not depend 
on e and reduces purely to phase rotations: 



Qs 



A 

1\A\ 



(6.5) 



where A can be obtained from Eq. ( 2.33 ). Therefore we can perform the trace in Eq. (6.4) over energy variables and 
the Keldysh matrix space, using Eq. fl6.5| ) and the formula |22f| 



dEtr [7V - -f l A (E + A {E - uj/2)] = 4w(iTj 1 ) 1 



(6.6) 



where the matrix II^ 1 is defined in Eq. (4.21). As a result, we obtain the effective action as a functional of the order 
parameter: 



o iGA . 

S A = - — . ,0 tr 



A| 2 n / 2tt 



A T (- w )o,nj 1 A( w ) 



(6.7) 



On deriving this equ ation we added the constant term Q s = 1 under the trace in Eq. (6.4). Employing Eq. ( 2.33 ), 
one can rewrite (6/7) in terms of the phase variables 9 as 



S A [9i] = ^ J ^{iu, [(e^ co S 2 ). a (e-^ sin 



9%) u - (e- ldl cos0 2 )-Ue tei ahx9 2 ) u 



2wcoth— {e Wl sinfe)-^ - * 1 sin0 2 ) u 



(6.8) 



The expression ( J6.8[ ) for the action makes it possible to relate the coefficient Ga in Eq. (S.4) with the Andreev 
conductance of the N-S interface. For this purpose let us suppose that the superconducting island is biased at some 
voltage V(t). Then the phase of the island will rotate with the speed 2eV. In the Keldysh formalism this corresponds 
to the rotation of its classical component, 



dffi 
~dT 



= 2eV. 



(6.9) 



To find the Andreev current, / = 2e^, where n is the number of the Cooper pairs on the island, one may use the fact 
that n and 9 are canonically conjugated variables, and thus, / = 2iejg. Translating this into the Keldysh formalism, 
we have, similar to |B2]: 



(J(t)) = le 



5\nZ 



59 2 {t) 



SZ 



02=0 



S9 2 (t) 



(6.10) 



02=0 



where the last equation follows from the normalization condition Z = 1 in the absence of quantum sources. 
Transforming Eq. (6.8) with the help of Eq. (6.9), we get 

S A = -eG A [ V{t)9 2 (t)dt + o[9 2 }. 



Substituting this action into Eq. (3.1C) and taking the functional derivative, we obtain 

(1(f)) - e 2 G A V(t). 

Hence, the Andreev conductance of the N-S interface in conventional units is equal to 



(6.11) 



(6.12) 



& A 



-G A . 



(6.13) 



This completes the proof of the physical meaning of the term Sa, Eq. (6.4), in the action. An alternative proof is 
presented in Appendix A, where we calculate the voltage noise at the N-S barrier (related to the conductance a a due 
to the fluctuation-dissipation relation) . Now we turn to the derivation of such a term in the simplest geometry of N-S 
contact. 
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C. Rectangular N-S contact 



1. Semiclassical solution within effective action method 



The simplest example of N-S contact is shown in Fig. || Superconductor (S) is connected to a normal reservoir (R) 
via a thin film of dirty metal (N) of length L x and width L y . The relevant energy scale in the N region is determined 
by the Thouless energy Eth — D/L x . As long as temperature, voltage and frequency of measurement are all small 
enough, ma,x(eV,T,ui) <C Erh, the Andreev conductance is just a constant: a a = Oj-Rd, where Rd = c~ 1 L x /L v is 
the resistance of the N region Below we first will show how to get such a result (described by the term (|6~|)) 
in the action) within the effective action method, and later on we turn to its generalizations, taking into account 
finite-energy effects as well as effects due to the interaction in the Cooper channel and due to the zero-bias anomaly. 




FIG. 4. Rectangular N-S contract. Superconductor (S) is connected to a normal reservoir (R) via a dirty metal film (N) of 
length L x and width L y . 



Consider the boundary action S-y defined in Eq. (6J) in the case when = Qs, = Q- In the present 

geometry it is convenient to treat the boundary as a line between N and S films, then 7 becomes the tunneling 
conductance per unit length. To obtain the Andreev term Sa, we need to expand the functional integral for Z up to 
the second order over 5 7 and average over fluctuations of the normal-metal matrix Q, i. e. over diffusion and Cooperon 
modes: 

S A =\{S*) = -^IAj\ytrQ s Q{Q,y)\ ^ (6.14) 

= / " d ^ / " ^ < tr Qs<r z r z W(0, yi) ■ tr Q s a z r z W(0, y 2 )) (6.15) 

d * Jo Jo 

Proceeding from Eq. Q6.14 ) to Eq. ( |6.15j ) we switched to the "rotated" representation of the Q-matrices defined in 
Eq. (4.4) and used the representation of Q in terms of the generators W as defined in Eq. (4.1C). Note that we neglect 
here any possible fluctuations of the superconductive matrix field Qs- The next step is to expand the field W(x,y) 
over the eigenfunctions ip mn (x,y) of the diffusion equation: 

W(X,V) =J2 W ^rnn(x,y). (6.16) 

The choice of the eigenfunctions is determined by the boundary conditions, which are: vanishing of current at the 
edges y = 0, y = L y and at the N-S boundary x = (the last condition is an approximation valid due to smallness 
of GtRd), and vanishing of the Cooperon amplitude at the boundary with the normal reservoir. These conditions 
result in the following set of eigenfunctions: 

2 cos q m x cos k n y n ( l\ tt 

iPmn{x,y) = —f= ==== , 1m = -r\m+-\, k n = —n, m,n = 0, 1,... (6.17) 

V(l + n ,0)L x Ly \ * J 

After the integration over y in Eq. fl6T5|) , only zero mode n = survives. Using Eq. ( |4.15| ) for the contraction rule 
and the fact that {Qs,t z } = 0, we conclude that the diffusion pairing in ( 4.16 ) does not contribute, whereas the 
Cooperon one does, so we get 

_ iir(~/L y ) 2 f deide 2 \ - tr(Qg(ei, ei)a z r z Q s (t 2 , ei)a z T z - Qsjei, £2)85(62, ei))(£>g T 2 n + z(ei + e 2 )<r z ) 

A \^L x L y J (27T) 2 {Dq 2 m) 2 + (£1 + £2 )2 

in^Ly) 2 fdEduj^ tr[Qs(u)A(E-)Q s (-u)A(E + ) - Q s {u)Qs{-u)]{Dq 2 m + 2iEk {E+)) 

2^ m„2 ^2 z AF 2 . ( 6 - 18 ) 



WvL x L y J (2tt)2 ^ {DqlY+AE* 
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where E± = E ± ^. Using the property Q 2 S = 1 one can verify that the contribution of the term with iEAo(E + ) in 
the numerator behaves as uj 2 for small uj. For this reason it can be neglected compared to the u-proportional kernel 
in Eq. (6.7). However, we will see below that the term of this structure determines the amplitude of the Josephson 
proximity coupling. Then, the term with Dq 2 n Qs{uj)Qs{—uj), being integrated over u>, produces an unimportant 
constant that can be omitted. As a result, one has 



S A = 



iAlL y ) 2 [dEdu^ Dq 2 m tr[Q s (m)A(E_)Qs(-u)A(E + )} 



\QvL x L y 



(27T) 



2 / j 



(Dq 2 n ) 2 +4E 2 



(6.19) 



In the limit E <C E^h = D/L 2 the sum over m can be easily calculated: 



1 



^ Dql 



TT 2 D 



i 



\) 2 



2D' 



(6.20) 



The region of energies E which are relevant for the integral in E q. (|6 .18 ) is given by E ~ max(cj,T). Thus, as long 
as the condition max(w,T) <C E^h is fulfilled, we can, using Eq. (S.20), convert (6.18) into the foreseen expression of 
the form (O): 



16 g L y 



(6.21) 



where Gt = "fL v is the total normal-state conductance of the interface, according to Eq. (6.3). Now, comparing with 
Eq. (6.4), one finds for the dimensionless Andreev conductance 



q(o) _ G\ L x 
A 9 Ly' 



(6.22) 



The above relation coincides with the known Q relation a a = ct^-Rd, since in the present geometry Ru = a^ 1 L x j L y . 

At higher temperatures or frequencies the ^-dependence of the denominator in Eq. ( |3.1£ ) becomes important and 
the simple invariant representation (6.4) is not valid anymore. However, one can still present the Andreev term in 
the action as a slight modification of Eq. (6.7): 



S A 



dw 



A| 2 n / 2tt 



■tr 



The function Ga(&) is given by 



G\ L eS (L x , uj, T) 



(6.23) 



(6.24) 



where 



L e s{L x ,uj,T) 



2D 

L x uj 



dE 



tanh ■ 



E+ 
2T 



tanh ■ 



E_ 
2T 



E 



{Dq 2 n 



AE 2 ' 



(6.25) 



In the limit max(w, T) <C Eth the result L e ff = L x is uncovered, whereas in the opposite case, max(w, T) ^> Exh, one 
gets 



tanh ■ 



E+ 
2T 



tanh ■ 



E- 
2T 



^E 1 0.95^/Dj2T, 



if T< w; 
if uj < T, 



(6.26) 



The number 0.95 in the above equation is the approximate value for (1 — 2 _3 / 2 )7r _1 / 2 C(3/2). 

Representation ( |6.23] ) can also be used in order to derive an expression for the nonlinear Andreev current I a (V) . We 
present this derivation in Appendix B, the result is that dc subgap current Ia(Y) — VGa(2bV) where the function 
Ga{<+>) is defined in Eq. J6.24 ). This relation between frequency-dependent linear response and static nonlinear 
response is due to the fact that static voltage V applied to the superconductor leads to the oscillation of its order 
parameter |A|e with the frequency 2eV. 

To summarize this subsubsection, we emphasize that within the usual semiclassical approximation the term de- 
scribing Andreev conductance is seen as a result of Gaussian integration over noninteracting Cooperon modes in the 
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N metal. Below we will take into account Cooperon nonlinearities which are due to interactions in the Cooper and 
direct Coulomb channels. Basically where are two different kinds of the interaction effects: the first one is due to the 
presence of electron-electron interaction in the Cooper channel (which is itself renormalized in a way dependent on the 
degree of disorder, as shown by Finkelstein) , whereas the second one is of the same nature as the zero-bias anomaly 
in the usual tunneling conductance |7|. These two effects are intrinsically different, as the first one is determined by 
diffusion modes with low frequencies u> -C Dq 2 and does not depend on the long-range tail of the Coulomb potential, 
whereas the second one comes from relatively high-frequency fluctuations with ui 3> Dq 2 and does depend on the 
actual behaviour of the Coulomb potential at large space scales (and thus it depends on the sample geometry as 
well). Below we start from the study of the first ( "Finkelstein's" ) effect and later on will include the effect of the 
zero-bias anomaly. The latter effects may indeed be strongly suppressed if a good conductor is placed near a dirty 
metal film, so that the Coulomb interaction in the film becomes screened on a relatively short distance (cf. Eq. ( 4.25| ) 
and Refs. 0,|TJ). 



2. The effect of interactions in the Cooper channel 

In this subsection we will show that the Cooper-channel interaction in the normal metal leads to a logarithmic 
(or weak power law) dependence of the Andreev conductance on frequency. To take it into account, one has to sum 
all inclusion of the vertex 5a into the Coo peron propagator. Such a summation can be effectively expressed as a 
renormalization of the tunneling action (3.1) due to the presence of the Cooper repulsion A. As a result, the tunneling 
conductance 7 gets renormalized down to 7(C) which should then be substituted into Eq. (6.24) for Ga- 



Qs 



g+co/2 E'+®/2 

>< Q 

E-a>/2 E'-(D/2 



FIG. 5. Correction to the interface transparency 7 due to interaction in the Cooper channel. 



Logarithmic corrections to the boundary term 
given by Eq. (|4.17|) : 



fh originate from pairing with the Cooper nonlinearity vertex S\ 



AS 7 = i(S^S\). 

The diagram is shown in Fig. |^, it contains one fast Cooperon mode connecting 5 7 and 



(S^Sx) — 



iirVyA 
16 



Tr r Qs<J z t z W ■ Tr{tr x ,Q}w g [T„W]i 



(6.27) 



(6.28) 



where Qs = uQsu, according to the general rule (4.4). 

Consider first the spacial structure of Eq. ( 6.28 ). The first trace corresponding to S 7 is taken along the N-S interface 
r, whereas the second one corresponding to S\ is taken in the bulk of the N region at a distance not larger than 
^/D/fi* from the interface. That is, from the point of view of slow variables, both S 7 and S\ are taken at the same 
point on the boundary T. In virtue of the coordinate integration in S\, the correlator (W'W) must be taken at zero 
momentum, so that the energy E becomes the only fast variable involved. Furthermore, due to the presence of the 
commutator of W' with t z under the second trace in Eq. (6.28), only Cooperon modes give a nonzero contribution 
(cf. Eq. (4.11)). Then, from the first trace one concludes that the whole average ( |6.28 ) does not vanish only in the 
sub-gap region, f2 < A. In other words, this reflects a trivial fact that the metallic tunneling conductance is not 
renormalized by interaction with Cooperons. In the deep sub-gap limit, the superconductive Qs is given by Eq. ( |6.5| ), 
it is independent of the "center-of-mass" energy E and depends on u> through the rotating phase 9(t). 

Employing the Cooperon-related part of the contraction rule (4.15), one obtains 



{S 1 S\j 



7T 2 7A 



dr 



duodEdE' 1 



32 J r "~ J (2tt) 3 E 
where il* < E < f2, whereas ui, E' < 51*, and the matrices 



tr(AB + Ao z Bo z - At z Bt z - Aa z r z Ba z r z )a z 



A = u E Q s (uj)u E a z T z , 

B = U E {a x , Q E> /2,E' +u /2\ u eo- Z T Z 



(6.29) 



(6.30a) 
(6.30b) 
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are taken at the same point at the interface (to logarithmic accuracy, all fast energies are equal to E). 
Using the property {Qs, r z } = and the relation Aq = ua z u, we transform Eq. ( 6.29| ) as 



TT 2 jA 

16 



dr 



dujdEdE' 1 



(27T) 



E 



tv{Q s ,{a x ,Q}}A (E). 



(6.31) 



The matrix Ao (E) in Eq. ( 6.31 ) can be replaced by 2 a 

by parity after integration over E. Employing the relations Qs oc Ik and {o+, a x } 



tanh ^ , since contributions of its other components vanish 



1, one arrives at 



(SjSx) — — 



TT"f\ f2 



In jy~ Trr QsQ- 



(6.32) 



This term effectively renormalizes the interface transparency; substituting Eq. ( 6.32j ) into Eq. (3.27), we obtain the 
RG equation for 7: 



£?7 



^7 = -7A. 



(6.33) 



A few comments are in order concerning this equation. First of all, we remind that it is valid in the sub-gap region 
C > Ca, where Ca = hi 1/Ar; for higher energies 7 remains unaffected. In the derivation we assumed that both ui and 
T are small compared to the running cutoff f2, otherwise the RG should stop at the scale max(u;,T). 

Eq. (6.33) contains the running coupling constant A(C) which is determined by Eq. ( |5.31 ) with Ao < replaced by 



A„ > — X g , a Cooper-channel repulsion in the N metal, defined at the energy scale r x . Looking at Eq. ( |6.33 ) one 



can think that it is valid for any geometry of the system since it arises as a result of integration over energy. This 
is however not the case and the RG exists in 2D only. The point is that the coupling constant A(£) logarithmically 
depends on the scale in 2D case only (cf. section |V Cj ). Recently a method to treat T c suppression for dirty SC films 
in the 1D-2D crossover region was developed p3| , which does not employ the existence of logarithmic RG equations. 
Substituting A(C) from Eq. ( |5.29| ), we write down the solution of Eq. (6.33) with the boundary condition 7(Ca) = 7o : 



7(0 



70 



1+ ^tanhA 9 (C- Ca) j <-oshA„k - i 



(6.34) 



Now this equation can be used to find the Andreev conductance modified by the interaction in the Cooper channel. 
In this regard we note that all integrals in Eq. (6. IS) are not logarithmic. Consequently, with logarithmic accuracy, 
it is sufficient to use the semiclassical expression for Ga, but with the renormalized barrier transparency: 



G A 



7^(_n^) (o) 
To 



(6.35) 



where the semiclassical value of G^' is given by Eq. (6.24) and the low-energy RG cutoff fi = max(D/L x , w, eV, T, t± 1 ) 1 



with 70 being the electron decoherence time in the N metal (all these quantities are assumed to be much below A). 
Substituting 7(C) from Eq. (6.34) one obtains 



Ga = — 



G\ L e s(L x ,u),T) 
9 



4(fi/A) 



2A„ 



[(1 + A n /A g ) + (1 - A„/A 9 )(ft/A) 2 ^ 



|2 ' 



(6.36) 



where A g is defined in Eq. (5.30). This equation is valid provided that the problem is effectively two-dimensional, 
cf. disc ussion after Eq. ( 16.33 ). This condition is satisfied as long as < L y , where Lq = y/D/£l. According 
to Eq. ( 6.36 ), in the low- Q limit the Andreev conductance acquires an anomalous power-law suppression with an 
exponent 2A 9 . The total power-law exponent, describing growth of Ga(Q) with the f2 decrease, is equal to 



XA 




(6.37) 
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D. Small SC island in contact with 2D film 



Here we consider another example of a N-S contact, namely, a contact of a small superconductive island of size d 
with a thin dirty normal film of linear size L ^> d, shown in Fig. Ifj. The edge of the film is connected to a normal 
reservoir. The semiclassical result for the Andreev conductance in this geometry can be inferred from the general 
relation a A l — a\Ru 1 where the metal resistance Rd is now a "spreading" resistance which grows logarithmically with 
the film size: Rd — m ht- If temperature or frequency is larger than the Thouless energy scale D/L 2 , the shortest 

of the effective lengths L^ ^ — ^J^pjj should be used instead of L. Logarithmic growth of the Andreev conductance 

with the space scale implies that G a itself becomes a running coupling constant subject to the RG procedure. Below 
we will derive and solve the RG equation for the Andreev conductance in the presence of the Cooper interaction 
renormalized by the Coulomb repulsion (still we will not touch here the effect of the tunneling DOS suppression; it 
will be considered below, in the next subsection). 




FIG. 6. Small superconductive island (S) of size d connected to a reservoir (R) through a dirty normal film (N) of size L d. 
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FIG. 7. Correction to the Andreev term Sa 



boundary vertices 5 7 : 



Similarly to Eq. (3.14), the correction to the Andreev term ( |6.4D in the action comes from the averaging of two 

(6.38) 



At scales larger than its size, d, the SC island can be considered as a point object. Then the boundary term (3.1) can 
be written in a local form: S-y = Tr QsQ(O), with Gt — Aj, where A is the area of the contact between the 



island and the film. The relevant pairing in Eq. (3.38), 

7T 2 G| 



(St) = 1 



16 



(6.39) 



is shown diagrammatically in Fig. ^. Here all energies are coupled to the slow matrices U, and the internal Cooperon 
momentum is the only fast variable. Integrating out fast modes with the help of the contraction rule ( 4.16| ), one 
obtains 



-kG\ f de±de2 dq 1 



I61/ J (2tt) 4 Dq 2 
where (all fields are taken at the point r = 0) 



Tr(AB — A(j z t z B<j z t z ), 



(6.40) 
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A={UQ s U- 1 ) ei , e2 a z T z , 
B = (UQsU- 1 ) i ^ 1 a z T z . 



(6.41a) 
(6.41b) 



Only the first term under the trace in Eq. (6.40) is important. Evaluating the logarithmic integral over q, one finds 

(6.42) 



Eq. ( 6.42 ) is valid at space scales L larger than the size of the island, d, i. e. for energies < ujd = D/d 2 . It shows 
that the RG procedure generates the term in the action of the form 



S A = -G A Tr(Q s Q) 2 , 
16 



(6.43) 



which reduces to Eq. (OA) when all fast modes of Q are integrated out and Q is replaced by A. The running constant 
Ga obeys the following RG equation: 



d( 4irg ' 



Gp 



(6.44) 



valid for £ > Q — In 1/Wr. 
In the semiclassical limit, 

above-mentioned result 



7 

Gel 
A — 



Gt/A is constant and integrating Eq. ( 6.44 ) over £ = 2\n{L/l) one reveals the 

G 2 T • A 1 

g-^-ln^. Beyond semiclassics, Eq. (6.44) becomes nontrivial since fluctuations which 
renormalize the tunneling conductance 7(C) must be taken into account. Below we will assume that cj^ ~ A and 
ther eby n eglect small corrections of the order of A 9 |Ca — Cd\- Then, substituting ^(() from Eq. ( |6.34 ) and integrating 
Eq. ( 6.44 ) over £ between Q and In we obtain the result for G^(O): 



GU(Q) 



Gp 



1 - {Q/A) 2 ^ 



47R7 (A 9 + A„) + (A 9 - A„)(f>/A) 2A « 



A general feature of the expression ( |6l5| ) is that G A (ty approaches a constant value, 

The semiclassical expression for the Andreev conductance, G c A \fl) 
increase of the relevant space scale Lq 



(6.45) 



in the low-f2 limit. 



4-n-g A„+A s 

G 2 

ji- In ^ , predicts growth of Ga(Q) with the 



i/D/n, due to the growth of the region where electron and Andreev-reflected 



hole interfere constructively. This expression for G c A l (fl) follows from Eq. (6.45) in the limit A„ — > and In <g; ^fg. 



Our result (6.45) demonstrates, that (in the present geometry when the current flow is effectively 2-dimensional) the 



Cooper-channel repulsion acts as if it imposes an upper limit for the time duration of such a constructive interference. 
The same qualitat ive b ehaviour would be seen in the absence of Finkelstein's corrections as well: taking first the limit 
g ^ In 2 — in Eq. ( |3.45| ) , one would find Ga ~ G\ / A-Kg\ n as Q — > (cf. Ref . Q] where a similar expression was used 
for Ga', note however that a numerical mistake was made in [Q, which lead to the overestimation of Ga by a factor 
of 2). Below we will see that the effect of the Coulomb-blockade suppression of tunneling is even more drastic, as it 
leads to the decrease of the Andreev conductance at fi — > 0. 



E. Effect of the zero-bias anomaly 

1. General treatment 



In the preceding consideration of the Andreev conductivity we neglected the effect of high frequency (Dq 2 <§; co) 
Coulomb fluctuations. They are responsible for the zero-bias anomaly (ZBA) in the usual tunneling conductance 
(calculated originally by Altshuler and Aronov 0), which, in 2D, leads to a suppression of tunneling by the relative 
order of g~ x In 2 1/flr. Later it was argued by Finkelstein |12] and shown (within semiclassical approach) by Levitov 



and Shytov |41| that the first correction must be exponentiated to get a result valid in the low-frequency limit as well. 
Recently, Kamenev and Andreev |22j rederived this result microscopically using the Keldysh cr-model approach and 



separating the Coulomb phase K according to Eq. (3.13). 

After the gauge transformation (3.13), the boundary tunneling term (pTj) acquires the form 
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S 7 = ^ 7 Tr r e iK " T 'QWe 



iK 12 Tzq{2) 



(6.46) 



where Kyi = K^ — K^ is the Coulomb phase difference across the interface. In our study of the Andreev conduc- 
tance, we will assume that the superconductor is connected to a low-impedance external environment which keeps 
fixed its average electric potential, and neglect the Coulomb phase K$ in the superconductor for the reasons explained 
in the last paragraph of Sec. IV C| . Hence, the ZBA is determi ned b y the fluctuations of the normal-metal phase K. 

On singling out the Coulomb phase by the transformation ( 3.13 ), the calculation of the ZBA becomes very sim- 
ple |22|: one has just to average the phase factors e in the boundary term (6.46) over the Gaussian fluctuations of 



K with the correlator (4.31). Indeed, (KiKj), being integrated over frequency and momentum, yields the aforemen- 



tioned g^ln 1/Qt. On the other hand, since the bulk action (4.5) depends only on derivatives of K, any {KiKj) 
pairing in the bulk contains an additional power of frequency or momentum, its contribution is less singular than 
In 2 1/VLt in the limit O — > 0, and therefore can be neglected. 




/ ✓ v \ 

FIG. 8. An example of factorization of the ZBA-type fluctuations for the Andreev conductance. The dashed zigzag line 



denotes the correlator {(e 



From this general observation it follows that in calculating the Andreev conductance the effect of the interaction in 
the Cooper channel is factorized from the ZBA effect. In other words, the full Cooperon which determines the subgap 
conductance is a product (in time domain) of the Cooperon without the ZBA (which was studied in subsections VI C 
and VI D| ) and the ZBA factor Z 2 (t) (cf. Eq. (7.15) below). Such a factorization is shown diagrammatically in Fig. |S|. 
Our previous analysis implied Z(t) = 1, i. e. weak ZBA at the scale ft ~ - r - 1 e 27r v / 9 relevant for the saturation of the 
RG flow for the Cooper repulsion constant A, cf. Eq. (5.31). We will use this factorization property explicitly in the 
study of the Josephson proximity coupling in Sec. VII. In the calculation of the Andreev conductivity we will rather 
follow another approach: the effect of the ZBA will be taken into account within the RG scheme together with the 
Cooper channel renormalization described above. 



2. Rectangular N-S contact 



High-frequency ZBA-type fluctuations of the Coulomb phase K in Eq. ( 6.46| ) lead to an additional suppression of the 
effective transparency 7 of the interface. This effect is different in the metallic (A < fi < r _1 ) and superconducting 
(f2 < A) regions. In the metallic case, the ZBA is governed by a single-electron tunneling, which is the only source 
of the transparency suppression, whereas in the superconducting case, the reduction of 7 becomes more pronounced 
due to a coherent tunneling of Cooper pairs, carrying the double charge. 



The renormalized value of 7 can be extracted from the boundary action ( 6.46 ) after averaging over fast fluctuations 
of the Coulomb phase K' in the metal: 



<5 7 ) = ^7Tr r ( e -^( tl )^Q 5 (t 1 ,i 2 )e^( t2 ^)Q(i 2 ,t 1 ), 



(6.47) 



where we neglected the Coulomb phase K$ in the superconductor as discussed above and used the fact that only the 
correlator of classical components of the phase K may be retained with logarithmic accuracy p2| ] . On eliminating fast 
degrees of freedom, the term S 7 will reproduce itself with the modified value of the transparency 7 (cf . |l0| ) . Further 
transformation of Eq. (6.47) depends on the structure of the matrix Qg in the Nambu and time spaces and on the 
geometry of the interface. 

In the metallic energy range, A < fi, the off-diagonal in the Nambu space components of the superconductive Qs 
can be neglected (cf. Eq. (3.11)), so that it commutes with the phase factor e lKl ^ Tz . Then the perturbative correction 
to 7 can be deduced from Eq. (6.47): 



27 



A7 = 7 



■ i<[K((tx)-iri(t2)] 2 ) 



(6.48) 



The fast phases K[ are correlated only at small times < fi^ 1 , while the time difference in Eq. ( 6.48| ) is large enough, 
I i 1 — £2 1 > due to an implicit e-dependence of Qs at the metallic region. Therefore, the cross average in Eq. ( 6. 48] ) 
can be neglected, and expanding the exponent one arrives at 



A 7 = - 7 ([^(0)] 5 



(6.49) 



The value of the correlator ([^(O)] 2 ) is not universal and depends on the setup considered. For the tunneling 
into the whole 2D plane (that corresponds to the geometry of a small SC island shown in Fig. |^) one can employ 
Eq. (pttl) : 



[Km 2 



1 plane 



a, 



dq 



Im 



V (q) 



IT J (2tt) 2 (Dq 2 - iuj)(2vV Q {q)Dq 2 + Dq 2 - iu)' 



(6.50) 



The leading double-log contribution to this expression comes from the region Dq 2 <C u> <C h'Vo(q)Dq 2 where the 
integrand is given by \/{gujq 2 ). Note that in this limit of strong Coulomb interaction, specific form of the potential 
Vo (q) drops from the integrand and enters the result only through the cutoff of logarithmic q- integration. For the case 
of the bare 2D Coulomb potential ( 4.24 ), one obtains 



plane — 




/27TfT 



8Tr 2 g 



In 



2 



-In 



2 ^0 



(6.51) 



where ujq = (2na) 2 /D. 

In the rectangular geometry of Fig 
the half-space was considered in Ref. 



a charge tunnels into the edge of the half-plane. A similar problem for 
Though a complete treatment of such problems is involved, a double-log 
asymptotics can be easily derived. To find it, one should use the fact that in the relevant region, Dq 2 <C u> <C 2i:aq, 
the screened Coulomb interaction is determined solely by the inverse density-density correlator fl45| whereas the bare 
Coulomb potential drops from equations as we have seen above. The density-density correlator can be obtained with 
the help of the eigenfunctions of the diffusion e quat ion with the proper boundary conditions. For the rectangular 
geo metry such eigenfunctions are given by Eq. ( 6.17 ). Then the value of ([if((0)] 2 ) can be obtained analogously to 
Eq. ( 6.51 ); it depends on the distance r from the edge of the half-plane and is given by 



mm? 



'half-plane 



1 



* 9 Jn, 



duj 



/2na 



dq 

2nq 2 



cos 2 (q x r) 



(6.52) 



where the infinite half-plane (L x , L y — > 00) is implied. For large r, the cosine squared in Eq. (6.52) can be substituted 
by its average value , 1/2, that leads to the infinite-plane result ( 6.51 ). For the tunneling directly into the edge (d = r) 
the correlator ( 6.52 ) appears to be two times larger. 

Hence the RG equation for the transparency 7 of a flat N-S boundary in the metallic energy range (C < Ca) reads 
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C-Co . 

2^ 2 ff 



•7> 



(6.53) 



where Co is defined at the scale luq, Co = lnl/woT. Being integrated over C from to In it yields exponential 
reduction of the tunneling conductance (which is two times stronger than the result jl2l|4 22 for the infinite plane): 



7(0 = 7o^m(C), 



where the factor Zm(() is given by 



Zm (C) = cxp - 



(C -Co) 2 -Co 2 

47r 2 g 



exp 



1,1, <r 



(6.54) 



(6.55) 



When the RG cutoff f2 becomes smaller than A, situation changes. First of all, Qs acquires an e-independent form 
( |6.5[ ), i. e. it becomes a (5-function in time domain, t\ =t2- It also anticommutes with e lKl ^ T ^ due to its structure in 
the Nambu space. In this limit, the Coulomb phase K can be considered as an additive correction to the SC phase 
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9 on the island, with the total effective phase # c ff = — 2K. Then the effect of the ZBA in the sub-gap limit can 
be attributed to the high-frequency fluctuations of 9 e g(t), destroying the Cooper pair coherence at large time scales. 
Repeating the steps that lead to Eqs. ( 6.48 ) and ( |5.49| ), one has instead 



A7 = 7 



3 -i([2^(t)] 2 ) _ 



= -2 1 ([K[(0)f 



(6.56) 



The effect of the ZBA anomaly should be taken i nto a ccount together with the renormalization of the transparency 
due to interactions in the Cooper channel, cf. Eq. (|6.33|) . Thus, we get the following RG equation for 7 valid in the 
sub-gap limit C > Ca: 



-A 7 - 



C-Co 



2 ~^ 
7T 9 



Solving this equation for C > Ca , we obtain a modification of Eq. (|6.34 ) 



7(0 



1 + tanh \ g (C - Ca) ) < '< >> ; I > A,, k - Ca ) 



(6.57) 



(6.58) 



where 7a = 7(Ca) = 7o^m(Ca) is the transparency renormalized by the ZBA in the metallic energy region A < Q < 
t _1 , and the multiplicative factor Zs(() accounts for the sub-gap ZBA effect: 



Z S (C) = cxp 



(C-Co) 2 -(Ca-Cq)^ 

2**g 



exp 



1 ] A 

2vr 2 5 n 



In- 



(6.59) 



So far we considered an infinite half-plane with L X1 L y — > 00. In analogy with the treatment in section VI C 2| , it 
might be tempting to substitute the renormalized transparency 7(ln^:) given by Eq. ( 3.58| ) into Eq. (6.24) to get 
the Andreev conductance Ga{Q) in the rectangular geometry. This however would be wrong. The point is that the 
expressio n (|6.52 ) for ([AT'(O)] 2 ) contains two coupled integrations over uj and q, with two frequency-dependent length 
scales, y 1 D /uj and 2iro~ /u>. The former is the usual diffusive length that shows how far a charge spreads during time 
u) . The latter is associated with the electric field propagation in a conducting medium. One logarithm in the 
double- log ZBA expression comes from the spacial region y/D/uj < R < 2na/oj. However for small enough ui, the 
length 2-Kaju! bec omes larger than the system size L x (it is assumed that L x < L y ). This finite-size effect can be 



accounted in Eq. ( 6. 52] ) by substituting max(w/27rer, 1/L X ) as a lower limit of g-integration. Physically this procedure 
means that we neglect external impedance of the circuit connected to our dirty film in comparison with the effective 
"spreading resistance" of the film (given by 7rr( Q, w here f unctio n r(C) is determined below in Eq. (|6.62|) ). It is then 
straightforward to generalize the RG equations ( |6.53| ) and ( 3.57 ) for finite systems. In the metallic energy region the 
transparency obeys 



while in the sub-gap limit one has 



|U-A 7 -2r(C) 7 . 



The function r(C) is defined as 



r(C) 




if il > ui a \ 
\iui a >VL> E Th ; 

if E Th > n, 



(6.60) 



(6.61) 



(6.62) 



where u> a = 2-ko~/L x and Ct/i = hil/i^T = 21n(i a //)- 

The solutions of Eqs. ( |6~60| ), ( |g7gT| ) can also be represented in the form ( |6.54| ), ( |6 . 58| ) with the modified functions 
Zm(C) an d Zs(C)- We will not present here the complete list of formulae for arbitrary values of uj a / A and w CT r but will 
focus instead on Zs{Q in the experimentally relevant case. If L x is measured in ^m then the energy ui a (in Kelvins) 
is given by 100g/L x , that appears to be greater than A for a reasonable experimental setup. Hence, in studying the 
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sub-gap frequency region, the lower g-cutoff is given by the inverse system size, L 1 , and any information about the 
Coulomb potential Vq (q) drops from the resulting expression. Solving then Eq. ( |6.61 ) one obtains 



Z s (0 =exp 



(C"Ca)[Ct/ 1 -(C + Ca)/2] 



exp 



1 A Afl 
tt- In — • In —75— 

2n 2 g Q E 2 h 



(6.63) 



Note that Q = max(i?T( 1 ,w, eV, T,t, ) and cannot be smaller then E 



Th- 



Using Eq. (6.63) and substituting 7(C) from Eq. (6.58) into Eq. (6.24), we obtain the resulting expression 



G. 



4(0/A) 



2A„ 



9 



L,< 



[(1 + x n /x g ) + (1 - \ n /x g )(n/A)^} 



exp 



1 



ir 2 g 



A Af2 

ln n- l »Ef h 



(6.64) 



where Gta = GtZm(C a) sta nds for the single-particle tunneling conductance at the scale A renormalized by the 
normal-metal ZBA. Eq. (6.64) is one of the main results of this paper, it shows that the original growth of Ga(^) oc 
Vl~ XA with the fl decrease (the exponent xa is defined in ( |6 .37 )) stops at ln(A/£!) ~ y^g, and at lower £1 the 
Andrccv conductance decreases due to the zero-bias anomaly. In the intermediate frequency range ln(A/f2) ~ ^fg 
both the Finkelstcin's effect and the ZBA effect are of the same importance, whereas in the infrared limit the ZBA 
effect is the most important one. The influence of the last effect upon the Andreev conductance was predicted by 
Huck, Hekking and Kramer |28| on phenomenological grounds. They considered N-S junction coupled to the model 
dissipative environment characterized by some impedance 2(u). We provide here a microscopic calculation of this 
effective impedance: Z(oj) 



7rr(ln ^r), with the function r(£) defined in Eqs. (6.62) and (6.68) below. 

*2 



Until now we considered the case of unscreened bare Coulomb potential Vo{q) = ^y-- In the presence of an 
additional (clean) metal gate, the Coulomb potential in the dirty metal layer changes according to Eq. ( 1.25 ). As a 
result, V r scr (0) = 47re 2 &, that modifies the law of propagation of the electric field, which now becomes diffusive with 
the effective diffusion coefficient = 8ni'e 2 b ■ D. The ratio D*/D = 2nib ^> 1, where K2 = Airi/e 2 is the inverse 2D 



screening radius. In a finite-size system, Eq. ( 6.62 ) should be modified as 



r (C) 




if §-n > 



e. 



Th, 



if f] > E Th > -g-VL; 

if E Th > n. 



(6.65) 



In an effectively infinite system (for Q > jy-Exh, when the propagating field does not have enough time to reach the 
edge of the s ystem), one of the two logarithms entering Eq. ( 3.51 ) becomes w-independent, cf. jyj. Then Ga is given 
by Eq. ( |6.64 ), provided that the last exponential factor is replaced by the power-law factor 



sen 2 

S \ 



(6.66) 



where 



— ln(87iW6). 



(6.67) 



3. SC island 



Now let us consider ZBA effects in the small island geometry shown in Fig. ^|. Let us first study the renormalization 
of the junction transparency 7. Depending on the relation between momentum and the inverse island size d , the 
geometry of the interface can be either flat, for qd ^ 1, or point-like, for qd <C 1. Then the renormalization of 7 
can be derived from Eqs. ( |6.49 ), (6.56), where the correlator ([i4T'(0)] 2 ) is given by Eq. (6.52), with cos 2 (q x r) being 



formally replaced by [1 + 9(qd — l)]/2. The resulting RG equations can be written in the form ( 3.6C| ), ( |6.61[ ) with the 
function r(£) defined as 



K0 = 



47r 2 g 



x < 



f 2(C — Co), ifft>27Rj/d; 

(,Th ~ Co, if 2ncr/d > D, > max(w (T , uj d ); 

C-Co, if u d > tt> u)<j\ 

C,Th - Crf) if "a > ^ > u>d; 

Crft - C, if min(w (T , u d ) > fl > E Th ; 

0, if E Th > n. 



(6.68) 
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The RG equations for 7(C) can be easily solved. Again, we will not present here a general solution depending on 
the relations between various energy scales, ETh{L) = D/L 2 , oj a = 2it<j/L, Ud, ft and 2ira/d. Rather we concentrate 
on the solution in the superconducting region, ft < A, assuming that u>d ~ A < w„. It can be written in the form 
( |6.58 ) with Zs(C) given by 



Z s {0 = exp 



(C-Ca)[Ct/,-(C + Ca)/2] 

47T 2 # 



exp 



1 1 A 1 
- —5- In — • In 

An 2 g ft 



An 



(6.69) 



The difference by 2 compared to Eq. ( 6.63 ) accounts for the difference between spreading over the whole plane and 

the half-plane. 

According to section VI D , in the island geometry, Ga entering the action as a parameter in Eq. ( 6.43| ) gets 
renormalized after eliminating fast degrees of freedom. In the absence of the ZBA, its renormalization comes from 
the averaging of two vertices 5 7 , cf. Eq. (6.38). With the ZBA effect taken into account, the Andreev action Sa will 
renormalize itself similar to the term S 7 in section VIE2j . The corresponding expression reads 



(S A ) = ^G A Tr(Q s (i 1 )e 2 ^( tl )^Q(i 1 ,t 2 )Q s (t 2 )e 2 ^^)^Q(t 2 ,i 1 )), (6.70) 
16 

that is written for ft < A when Qs given by Eq. ( |6.5| ) is local in time. On averaging over fast K[, one finds for the 
correction to Ga- 



AG a = G A [ e -i<[2^(ti)-2if;(*2)] 2 > _ 

Expanding the exponent and omitting the cross averages, one arrives at 

AG A = -±G A ([K' 1 {<S)] 2 ). 



As a result, we obtain the following modification of the RG equation (6.44): 

dG A A 2 j 2 



4r(C)G, 



(6.71) 



(6.72) 



(6.73) 



Taking 7(C) from Eq. (6 
conductance: 

G A (0) = 



d( ing 

and integrating this differential equation, one obtains the solution for the Andreev 



r 2 



1- (ft/A) 



2A„ 



4ng (A 9 + A„) + (A S -A„)(ft/A) 2 ^ 



exp 



1 A AO 



(6.74) 



where Gta is the normal-state tunneling conductance at the scale A renormalized by the ZBA. The coefficient in 
front of the double logarithm in the ZBA exponent is four times larger than the one for the single electron tunneling, 
due to the doubled charge of a Cooper pair. Expression (6.74) again shows that the ZBA effect upon Ga is described 
by the separate factor Z s (() as it should be due to the factorization property discussed above. With the frequency 
decrease, G^(ft) first grows logarithmically and then decreases as in the case of the r ectangular N-S contact. In the 
presence of a screening metal gate, the ZBA factor [Z s "} 2 = {n/A) x ^ 2 (cf. Eq. (B)). 



VII. JOSEPHSON PROXIMITY COUPLING 



A. General treatment and the RG equation 



In this section we first rederive semiclassical expressions for the proximity coupling between superconductors, 
separated by dirty normal metal, and then generalize them, taking into account quantum fluctuations, similar to the 
way it was done above for the Andreev conductance. The term in the effective action, which is responsible for this 
coupling, can be written in the form 



1 



(7.1) 



where superscripts ^' and ^ refer to two superconductive banks or islands. Using the low-energy representation (|6.5|) 



for Qs, and neglecting the Coulomb phase factors enp(iK(t)T z ), one rewrites Eq. (7.1) as 
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Sj = -IE j J sin 01 (*) sin e 2 (*)dt, (7.2) 

with 0i and 02 being the classical and quantum components of the phase difference between the superconductors, 
= — 6>( 2 '. The meaning of the term (7.2) becomes transparent in the initial basis before the Keldysh rotation 
( [2.30 ): it yields the usual expression for the Josephson coupling energy, — Ej cos©. Employing Eq. ( 6.10| ) one obtains 



the standard relation between the Josephson current and (the classical component of) the phase difference (with 
dimensional units restored): 



2eE 



- sin 0i. 



(7.3) 



The term (|t|) can be derived in the way very similar to the derivation of the expression ( J6.4| ) in section |VI H : 
one should consider the cross term in the perturbative correction to the action, 



(i) 



where Si?' is the 



boundary action ( |6.1[ ) at the interface with the superconductor (j), and the average is taken over fluctuations of the 
normal- metal matrix Q. 
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FIG. 9. Two superconductive islands (S) of size d with separation i? S> d coupled via a dirty normal film (N) of size L. 




FIG. 10. Two superconductive terminals (S) coupled via a dirty normal film (N) of sizes L x and L y 



The Josephson coupling energy Ej depends on the geometry of the system, and below we will find Ej for two 
geometries shown in Figs. ^ and |Io|, that are natural counterparts of Figs. ^ and |J discuss ed a b ove w ith respect to 
the Andreev conductance. In both cases Ej is given by the Fourier transform (cf. Eqs. (7.6), ( 7.11 )) of the zero- 
frequency Cooperon amplitude. The latter, C(q,ui) = J7(C), logarithmically depends on Q = max(Dq 2 ,uj,T). Within 
the semiclassical approximation B it is given by an integral over the "center-of-mass" Cooperon energy: 



4^7 



[F(E+) + F{E-)]EdE 
(Dq 2 ) 2 + 4E 2 



8z/ n 



(7.4) 



To t ake into account the effect of interaction in the Cooper channel that makes Gt = A~f scale-dependent, we replace 
Eq. (7A) by the following differential RG equation: 



dJ(0 , 2 7 2 (0 



(7.5) 



This RG equation coincides (up to a numerical factor) with the RG equation (6.44) for G a in the small island geometry. 
In the study of the Andreev conductance, such an RG equation appeared as a result of the ^-integration and thereby 
hold only for the geometry of Fig. ^. In the case of the Josephson coupling, Eq. ( ^5] ) comes from the ^-integration 
and thus is insensitive to a particular geometry holding for both setups shown in Figs. ||, [l0[ 
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Although the RG equations ( |6.44 ) and (7.5) look similar, the answers for Ga and Ej are quite different. The main 
difference is that J is not a final quantity, it should be Fourier transformed to get the Josephson energy Ej. As a 
result, only the lowest spacial modes effectively contribute to Ej, which appears to be exponentially suppressed ||] 
due to loss of coherence if temperature or frequency of phase fluctuations exceed the Thouless energy. Note, that in 
the case shown in Fig. ||, the relevant Thouless energy scale is determined by the inter-island distance R; we denote 
this energy as ui a ; = D/R 2 to discern it from Eth(L) = D/L 2 which will also enter results for the ZBA factor for the 
two-island setup. Below we will assume that both T and u) are much lower than Exh for the rectangular geometry, 
and much lower than lor for the island geometry. 

Eq. ( |7.5D is sufficient to find the proximity coupling energy in the absence of the zero-bias anomaly effects (which 
will be incorporated in the next subsection). We start from an example of two small SC islands of size d, separated 
by the distance R^> d (see Fig. pf). In this case 



(7.6) 



where C, q = hi(l/ Dq 2 r), and the use of the continuous Fourier transform implies that the distance between islands is 
much shorter that the total size of the film, L~^> R. Using the identity Jo (a;) = ~ -§^[xJi{x)}, and integrating by parts 
we obtain 



1 f°° 

Ej(R) = — / dqJ^qR) 
kR Jq 



Since J{C, q ) is a very slow function of q (provided that g ^> 1), the ^-integration converges rapidly near q 
Then we may take J at C, q = In and perform the remaining integration that yields 



Ej{R) 



1 dJ{Q 



ttR 2 d( 



C=ln # 



Combining Eqs. ( |7.q ) and (7.5) with Eq. ( 6.34 ), we obtain 



Ej{R) 



8nuR 2 [(i + \ n /\ g ) + (1 - A„/A g )(a/i?) 4A » 



(7.7) 
R-\ 

(7.8) 
(7.9) 



Thus, the effect of repulsive interaction in the Cooper channel is to produce an anomalous power-law exponent 
Xj = 2 + 4A 9 , describing decay of the Josephson coupling at long distances: Ej ~ R~ XJ . 

Now we proceed to the case of the rectangul ar ge ometry shown in Fig. [l0|. In this case the set of eigenfunctions 
with the proper boundary conditions reads (cf. ( |6.17| )) 



ip mn (x,y) 



2 cos q m x cos k n y 



+ <5 m ,o)(l + 3n,o)L x L y 



L.r 



m, 



m, n 



0,1,. 



(7.10) 



The difference with the case of two islands is that now one has to sum over the single component q m of the wavevector 
only (cf. section VIC): 



Ej{L Xl L y ) — 



oo 



y m— — oo 



(-l) m C(q m ). 



(7.11) 



Eq. ( |7.5| ) for J = const • Ga has been already solved in section VI D ; substituting the solution into Eq. ( 7.11 ) and 
using the equality X)m=-oo(~ 1)'™ = we transform Ej to the form 



Ej(L x ,L y ) 



G 2 ^ 



E 



(-ir 



L x L y \ g - A„ ^ (X g + A„ ) + (X g - A„ ) ( E T h tr 2 / A ) 2 A s 



(7.12) 



where Etk = D/L 2 . Consider this sum and write it as X)m=-oo( — -*-)"* fm- F° r A g <C 1, f m is a very slow function for 
m 7^ but has a cusp at m = 0. Then the sum can be well estimated as /o — /i and is equal to 



(X g -X n )(E Th /A) 2 ^ 



(X g + X n )[(X g + An) + (A s - X n )(E Th /A) 2 ^} 



(7.13) 
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As a result, we obtain 



Ej(L x , L y ) — 



G 2 ^ 



4A„ 



8is{\ g + Xn)L x L y (1 + X n /Xg) + (1 - X n /X g )(U/Lx) 



4A„ 



(7.14) 



In the above expression Gt is the total normal-state tunneling conductance proportional to the width L y of the 
contact. In the rectangular geometry a non-trivial exponen t x'j = 1 + 4A g enters the Ej(L x ) dependence: Ej oc L x Xj . 



In the case of weak disorder, A g ln(L x /^) <C 1, the result (7.14) reduces to the one obtained in pj, whereas Eq. ( |7.9| ) 
reduces to the form used in [Q (an extra factor 4/tt 2 in the Eq. (18) of [Q is due to the difference in the definition 
of the dimensionless tunneling conductance Gt)- Both of the above results (|f~~ 
zero-bias anomaly effects. Now we proceed to take them into account. 



(7.14) were obtained neglecting the 



B. Proximity coupling in the presence of the zero-bias anomaly 



As it was explained in Sec. VI E 1 , the effect of the Coulomb phase K(t) fluctuations (which are due to the long-range 



Coulomb potential) can be separated (factorized) from the effect of the interaction in the Cooper channel. Below we 
will use this factorization explicitly. 

Consider first the two island geometry of Fig. M. Using the gauge-transformed form ( 3.46 ) of the boundary action 
term, and repeating the steps that lead to Eqs. ( |7.5|) and (7.6), we obtain now the Josephson coupling energy Ej(R) 
as a time-domain convolution of the Cooperon and the ZBA exponent: 



Ej{R) = / dtZ 2 j{t,R) 



dq 



jqR 



C(q,t), 



(7.15) 



where C(q,t) = j ^J{C,)e lu)t is the time-domain Cooperon amplitude, and £ = lnl/Qr with £1 = max(Dq 2 ,Lu). 
The ZBA factor Zj(t,R) is defined as (cf. Eq. flS.SlD) 



lnZj(t,R) = ([^(t, (0,0)] 2 ) 



1 



o/D 



dq 



n 9 JO W Jma,x(uj/2-KtT,l/L) 1 



1 — cos cut J n (qR)) 



(7.16) 



where L ^> R is the size of the metal film. Here we assume that u>d ~ A, and include all ZBA fluctuations with 
frequency uo A into redefinition of the normal-state tunneling conductance Gta at the scale A. In the absence of 
the ZBA, Zj = 1 and Eq. Q7.15| ) reduces to Eq. (ffT^), with C{q) = C(q,u = 0). The lower limit of the g-integration 
depends on the relation between the field spreading distance, 2tt(t/uj, and the sheet size, 1/L. To be more specific, 
we will assume, following section VIE, that A/27rer -c L. Then the ZBA exponent is equal to 



Zj(t, R)=exp[--^- In — ■ In A 



E Th (L) 



(7.17) 



where ujr = D/R 2 , ExhiL) = D/L 2 , and it is assumed that t ~ u^ 1 . Now substituting J{Q from Eq. ( [7.51 ) into 
Eq. d7.15Q , using the trick that led to Eq. §T§ and the identity dJ{Q/dC, q = 6{Dq 2 - \uj\)dj(( q )/d( q , and performing 
a trivial integration over ui, we obtain 



Ej(R) 



TT 2 R 



dtZ 2 j(t,R) / dqJi(qR) 



dj(( q ) sin DgH 

d( q t ■ 



(7.18) 



Both momentum and time integrals converge fairly well in the vicinity of q 
result, to logarithmic accuracy we get 



Ej{R) 



1 



R 1 and t ~ ujp} respectively. As a 



(7.19) 



irR 2 d( 

Note that the same result could be obtained within the RG approach with the help of the equation similar to Eq. ( 3.73 ) 

A 2 7 2 



dj 



(7.20) 
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with r(C) = -(1/2) din Zj/d(. 

Finally, we obtain the Josephson coupling energy for two islands: 



Ej{R) 



r 2 



4A„ 



8nvR* [(i + X n /X g ) + (1 - X n /X g )(U/R) 



4A„ 



• exp 



—5— In — • In - — — 

K z g £a £aR 



(7.21) 



The presence of the total size L of the film in the above expression for the inter-island proximity coupling may appear 
to be unexpected. It originates from the fact that electric field propagates much faster than electron density: during 
diffusion time R 2 /D corresponding to the inter-island distance, the electric field propagates the distance L(R) = K2R 2 , 
where K2 — 47we 2 is the inverse 2D screening length. The effective electric impedance that determines the ZBA factor 
does depend on the distance L(R) as long as it is shorter than the total size L. If the dirty film is of metallic origin, 
k^ 1 is of the order of Angstroms, so L(R) easily exceeds any reasonable size of the film even for R in submicron 
range. In this case the effective impedance is determined by the film size L. All these considerations do not hold for 
the case of screened Coulomb potential ( 4.25| ); in that case the exponential factor in Eq. (7.21) should be replaced by 
(U/R) Xz , cf. the end of Section VI. 

Consider now the ZBA effect for the rectangular geometry (see Fig. nffl) . In this case one has 



Ej(L x ,L v ) = -J— f dtZ 2 (t,L x ) V (-l) m C(q m ,t), 

L x L,y J _ 



(7.22) 



m— — 00 



cf. Eqs. ( [fill ) and ( ffl^) . A gain, the relevant scale for the i-integration is given by E~l = R 2 /D. Then one can 
integrate over t keeping Zj(t) = Zj(E^) constant that amounts to multiplying the previous result J7TI41 ) by Z 2 (£"*), 
provided that Gt is substituted by Gta- The ZBA factor for a flat N-S interface is given by 



Z j( E Th' L x) = ex P ( - 



2n 2 g 



E 



Th 



where Eth — D/L 2 . Therefore we obtain 
Ej(L X} Ly 



r 2 



2(U/L X ) 



4A„ 



8^(A 9 + X n )L x L y (1 + X n /Xg) + (1 - X n /X g )(U/L x ) 4X * 



■ exp 



Xn 2 hi 



For the screened Coulomb potential (parallel metal gate), the last factor in (7.24) is replaced by (^a/^3 
finally, that both expressions (7.21) and (7.24) refer to the zero-temperature limit T <C hD/R 2 , hD/L 2 . 



(7.23) 

(7.24) 
) 2x *. Note, 



VIII. DISCUSSION 



We have developed field theory functional formalism of the Keldysh type for disordered superconductors. This 
approach provides a regular method for the treatment of all kinds of quantum fluctuations beyond the standard 
semiclassical theory of superconductivity []l8|-^of which can be understood as the saddle-point approximation of our 
theory The theory is formulated in terms of the local (in space) matrix order parameter whose components corresponds 
to the retarded/advanced and Keldysh Green functions. The main advantage of our approach with respect to the 
standard Matsubara replica formalism |]l2| is that it provides a possibility to treat nonequilibrium problems. 

General formulation of the theory involves 8x8 matrices Qt,t', composed as the direct product of the 2x2 blocks 
corresponding to the Nambu, time-reversal and Keldysh spaces. However, in order to present our approach in the most 
transparent form, we restricted here our specific calculations to the case of spin-independent interactions, and neglected 
usual weak-localization conductivity corrections, that makes it possible to use reduced 4x4 matrices. We focused 
here on the fluctuational effects specific for the dirty superconductive films, those relative magnitude is known [jll 12 
to be of the order of (1 / g) In 2 (L / 1) (in the two-dimensional limit), and thus can be considered independently of the 
normal-metal weak-localization and interaction corrections [^JQ| which are cx (1/g) ln(L/Z). There are two physically 
different types of these effects: i) Coulomb-induced repulsive contribution to the electron-electron interaction in 
the Cooper channel (coming from the "diffusive" frequency region lj -C Dq 2 ), and ii) reduction of the averaged 
single-particle density of states due to high-frequency (u> ^> Dq 2 ) fluctuations of electric potential. We have shown, 
following the approach developed in p2], that these two effects can be separated at the non-perturbative level, by 
"gauging away" phase factors induced by the fluctuating long-range electric potential. The first effect was treated 
by the renormalization group method within the single-loop approximation. Being applied to the case of uniformly 
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disordered superconductive films, this method essentially reproduces the well-known results by Finkelstcin (TT 13 for 
the Coulomb-induced suppression of the superconductive transition temperature T c . The second effect is formally 
similar to the "infrared catastrophe" of quantum electrodynamics; in our technique it is taken into account in an 
essentially exact way, similar to the normal- metal case treated in 22 1 . This second effect (analogous to the zero-bias 
anomaly of the tunneling conductance), contrary to the first one, does depend on the long-range part of the bare 
Coulomb potential, and, thus, can be modified by the presence of an additional external screening (provided, e. g., 
by a metallic gate nearby the dirty film). 

It was argued previously (cf. (ll|) that the fluctuations responsible for the zero-bias anomaly effect are specific for 
the single-particle density of states, due to gauge-noninvariance of the single-particle Green function, whereas the 
same fluctuations do not contribute to other (gauge-invariant) physical quantities like T c . We have calculated two 
important characteristics of mesoscopic S/N structures, the Andreev subgap conductance and the Josephson proximity 
coupling, and found that they are affected by the zero-biased anomaly effect as well as by the Finkelstein's effect. To 
be more specific, we have calculated linear Andreev conductance Ga(Q) in the presence of the tunneling barrier (N-I-S 
structure), for the two different geometries: i) small SC island "sitting" on a large-area dirty metal (Fig. ^), and ii) 
rectangular no rmal film b etwe en superconductive and normal reservoirs (Fig. [|). Our results for those two cases are 
given by Eqs. ( 6.74 ) and ( |6.64 ) correspondingly. In these formulae f2 stands for the infrared cutoff of renormalization 
procedure, i. e. Q = max(T, oj,t^,,D/L 2 ). A detailed generalization of these results for the steady-state nonlinear 
Andreev current Ia(V) can be obtained with the use of formulae from Appendix B; qualitatively the behaviour of 
G A (V) = dI A /dV can be found by substituting 2eV for ft in Eqs. ( |6.74| ) and fl6.64| ). 

The second new physical quantity we have studied is the Josephson proximity coupling Ej as a function of the size 
of the normal region between two SC contacts. It was calculated, in the low-temperature limit T <C Eth, for the 
case of two small islands in contact with dirty normal metal (Fig. ^|), as well as for the rectangular film between two 
SC banks (Fig. [ll]). The electron-electron repulsion in the Cooper channel leads to the appearance of an anomalous 
scaling for Ej. In particular, for two small SC islands Ej(R) cx i?~ 2 ~ 2 / 7r v / ff , whereas for the rectangular contact the 

Josephson energy decays with the film length as Ej(L x ) cx L x 1 2 '^^, I n addition, the "ZBA effect" adds to the 
above behaviour a log-normal suppression factor exp[— (A/n 2 g) hi 2 (L x /^)] for the rectangular contact, and a similar 
factor for the isl and g eometry. Full results for both geometries are given, at the low-temperature limit T — > 0, by 
Eqs. ( 7.21 ) and ( 7.24 ). Screening of the long-range Coulomb interaction in the film by an external gate makes the 
ZBA suppression weaker, it reduces then to the power-law fa ctor ( £a/ Lx) 2x * and ((,a/R) x * for the rectangular and 
island geometries correspondingly, where x z is defined in Eq. ( |3.67| ). 

In the present calculations we did not took into account the existence of a finite decoherence time in normal 
metals at T > 0. Qualitatively, the effect of nonzero r^ 1 is to suppress the central peak in I(V) predicted within the 
semiclassical theory, i. e. it acts in the way similar to the effect of quantum corrections we have studied in this paper 
(we emphasize that these quantum effects have nothing to do with "zero-temperature decoherence" ) . That is why the 
decoherence rate t^ 1 extracted from the data on the Andreev conductance might be overestimated unless quantum 
fluctuations are taken into account in the data analysis. Another unsolved theoretical problem is to go beyond the 
lowest-order expansion over the tunneling conductance Gt, that becomes necessary if Gt approaches or exceeds the 
conductance of the metal region. 
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APPENDIX A: ALTERNATIVE PROOF OF EQ. (6.13) 



In this Appendix we present an alternative proof of Eq. ( |6.13 ) with the help of the fluctuation-dissipation theorem. 
Expanding the action (fx8) to the second order in the phase variables, we obtain 
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Hence, in the Gaussian approximation, the correlator of phases is given by 

2(nj« 



(Al) 



(A2) 



Employing Eq. (6.9) relating the phase of the island with the applied voltage, we get for the spectral voltage correlator: 

(A3) 



that, according to the Nyquist formula, results in Eq. (3.13) 



APPENDIX B: NONLINEAR ANDREEV CONDUCTANCE 



To derive nonlinear Andreev current I A (V) we start from Eq. (6.23). Tracing over the Nambu space we reduce the 
expression for the action to the form similar to Eq. (S.8), with Ga being replaced by G A (u>). Calculating the current 
following Eq. (6.10), we find 



— iu>{ti — t) 



/ e ifli(ti)-ifli(t)^ _ ( e -*0i(ti)+-i0i(*)) _ 



2i coth — 



{e iei{tl) e 2 {t 1 )e~ tei{t) ) + (e- iei ^6 2 (ti)e m ^) 



(Bl) 



Angular brackets in Eq. (Bl) mean averaging over quantum dynamics of the phase 6{t). In the case when fluctuations 
of 9 can be neglected, the second line in Eq. (Bl) should be omitted, whereas in the first line we put e^ittv-Wiw = 
e 2ieV(*i-t)_ As a result, integration over ti produces 2ttS(uj ~~ 2eV), and one obtains 



I A = e 2 VG A (2eV). 



(B2) 



The function G(uj) in Eq. (B2) should be determined with the account of fluctuational corrections due to the interaction 
in the normal metal, as discuss ed in Sees. VIC, VI D . The effect of the zero-bias anomaly also can be taken into account 
along the same lines, cf. Sec. VIE . Ho wever, the RG method employed there does not determine the dependence 
of the ZBA factor Z S (Q) (cf. Eq . d^59| ) upon the ratio 2eV/T. To find this dependence, we use the factorization 
property discussed in Sec. |VI E 1 and note that the low-energy ZBA effect can be accounted for by the replacement 
9(t) — > c s{t) — 9(t) — 2K(t), cf. discussion after Eq. ( 6.55| ). After such a replacement the correlation functions 
entering Eq. (B4) can be expressed via the correlation matrix 



V{UJ) - I V A {u) 



(B3) 



where, similar to Eqs. (|4.3C| ) and ( fl.3l|) , the matrix elements T>ij(uj) = — 2i(*i(o;)*+ (—w)), and %>(t) = e r 2iK{t \ As 
a result, the whole expression for the current can be written as 



-^wGaH iI> K (2eF-w)+2coth — lmV R {2eV - w) 
2tt L 2T 



(B4) 



Note that iD K (w) is a real positive function. Below it will be convenient to rewrite Eq. (JB4|) via the forward/backward 
correlation functions U > ( < }(uj) defined as 



V K = V> 



V R -V A = T>> -T> < . 



(B5) 



In the equilibrium state, the detailed balance relation T> < (lo) = e~ UJ ^ T V > (u>) is valid for bosonic correlation functions, 
in addition to the general relations (B5). In terms of the original (before the Keldysh rotation ( 2.3C )) bosonic variables, 
T> > (t) = — z(\E f (t)\E f+ (0)). The correlators labeled by the index ">" contain fields taken at the forward branch of the 
Keldysh contour C (for the sake of brevity we will omit the corresponding index "1"). 
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If correction to Ga(w) from the interaction vertex A in the Cooper channel can be neglected, it can written as 



dE 

U! 



, E + i E - 

tanh tanh 

2T 2T 



C(E), 



where E± — E ± u/2, and 



f 



-2iE + Dq 2 



(B6) 



(B7) 



is the Cooperon amplitude. The sum in Eq. (B7) goes over Cooperon eigenmodes (cf. derivation in Sec. VIC 1). 
Combining Eqs. (B4)-flB7|), we obtain finally 



Ia(V) 



C{E)dE 



00 dEy_ 



, E-E'/2 + eV , E + E'/2~eV 
tanh — tanh — 



2T 



2T 



V>(E')- 



-2eV/T 



(E'-2eV)/T- 



(B8) 



The result (B8) is equivalent to the one obtained in |28|] by the method of analytic continuation. In the T — > limit 
expression (B8) simplifies and can be transformed to the following form for the differential conductance: 



dV 



2^2 r 2eV 



e 2 G, 



/zev 
dEV(E)C(eV - E/2), 
-OO 



(B9) 



where V(E) = (i/2n)V>(E). 

Note finally that the results (B8), (B9) are valid in the situation when the order parameter phase 9(t) is subject to 
quantum fluctuations as well; in this case the correlation function T> > is defined as r D > (t) = —i{e ie ^^{t)e~ l9 ^°'^ + {Q)), 
where 9{t) = 0(t) - 2eVt. 



A. F. Volkov, A. V. Zaitsev and T. M. Klapwijk, Physica C210, 21 (1993). 
Yu. V. Nazarov, Phys. Rev. Lett. 73, 1420 (1994). 

C. W. J. Beenakker, Rev. Mod. Phys. 69, 731 (1997). 

L. G. Aslamazov, A. I. Larkin and Yu. N. Ovchinnikov, ZhETF 55, 323 (1968) [Sov. Phys. JETP 28, 171 (1969)]. 

E. Abrahams, P. W. Anderson, D. C. Licciardello and T. V. Ramakrishnan, Phys. Rev. Lett. 42, 673 (1979). 

L. P. Gorkov, A. I. Larkin and D. E. Khmelnitsky, Pis'ma ZhETF 30, 248 (1979) [JETP Letters 30, 228 (1979)]. 

B. L. Altshuler and A. G. Aronov, in Electron- Electron Interactions in Disordered Conductors, ed. A. J. Efros and M. Pol- 
lack, Elsevier Science Publishers, North-Holland, 1985. 

B. L. Altshuler, A. G. Aronov and P. A. Lee, Phys. Rev. Lett. 44, 1288 (1980). 
S. Maekawa and H. Fukuyama, J. Phys. Soc. Jpn. 51, 1380 (1982). 

H. Takagi and Y. Kuroda, Solid State Comm. 41, 643 (1982). 

A. M. Finkelstein, Pis'ma ZhETF 45, 37 (1987) [JETP Letters 45, 46 (1987)]. 

A. M. Finkelstein, Electron Liquid in Disordered Conductors, Vol. 14 of Soviet Scientific Reviews, ed. I. M. Khalatnikov, 
Harwood Academic Publishers, London, 1990. 
A. M. Finkelstein, Physica B 197, 636 (1994). 

Yu. N. Ovchinnikov, ZhETF 64, 719 (1973) [Sov. Phys. JETP 37, 336 (1973)]. 

K. B. Efetov, A. I. Larkin and D. E. Khmelnitsky, ZhETF 79, 1120 (1980) [Sov. Phys. JETP 52, 568 (1980)]. 

K. B. Efetov, Adv. Phys. 32, 53 (1983); Super symmetry in Disorder and Chaos, Cambridge Univ. Press, New York (1997). 

L. V. Keldysh, ZhETF 47 1515 (1964) [Sov. Phys. JETP 20, 1018 (1965)]. 

A. I. Larkin and Yu. N. Ovchinnikov, in Nonequilibrium Superconductivity, eds. D. N. Langenberg and A. I. Larkin (Elsevier, 
New York, 1986). 

J. Rammer and H. Smith, Rev. Mod. Phys 58, 323 (1986). 

C. J. Lambert and R. Raimondi, J. Phys. Cond. Matt. 10, 901 (1998). 
K. Usadel, Phys. Rev. Lett. 2 5, 507 (1970). 



A. Kamenev and A. Andreev, cond-mat/9810191, submitted to P hys. Rev. B. 

C. Chamon, A. W. W. Ludwig and C. Nayak, |sond-mat/9810282j , submitted to Phys. Rev. B. 



A. G. Aronov, A. S. Ioselevich, Pis'ma ZhETF 41, 71 (1985) [JETP Letters 41, 84 (1985)]. 



38 



C. Chou, B. Su, B. L. Hao and L. Yu, Phys. Rep. 118, 1 (1985). 
S. Babichenko and A. N. Kozlov, Solid State Comm. 59, 39 (1986). 
L. Horbach and G. Schoen, Annalen der Physik 2, 51 (1993). 

Huck, F. W. J. Hekking and B. Kramer, Europhys. Lett. 41, 201 (1998); |sond-mat/9709344 . 

N. Popov, Functional integrals in quantum field theory and stastistical physics (Reidel, Dordrecht, 1983). 

P. Gorkov, ZhETF 34, 735 (1958) [Sov. Phys. JETP 7, 505 (1958)]. 

Nambu, Phys. Rev. 117, 648 (1960). 

Altland, B. D. Simons and D. Taras-Semchuk, Pis'ma ZhETF 67, 21 (1998) [JETP Lett. 67, 22 (1998)]; 



extended 



version: cond-mat/9807371 . 

G. Eilenberger, Z. Phys. B, 214, 195 (1968). 

A. A. Abrikosov, L. P. Gorkov and I. E. Dzyaloshinskii, Methods of Quantum Field Theory in Statistical Physics, Prentice 
Hall, New York (1963). 

B. N. Narozhny, I. L. Aleiner, B. L. Altshuler, [:ond-mat/9903239 . 

A. I. Larkin and A. B. Migdal, ZhETF 44, 1703 (1963) [Sov. Phys. JETP 17, 1146 (1963)]. 
I. L. Aleiner and B. L. Altshuler, unpublished. 

V. G. Vaks and A. I. Larkin, ZhETF 40, 1392 [Sov. Phys. JETP 13, 979 (1961)]. 

M. Yu. Kupriyanov and V. F. Lukichev, ZhETF 94, 139 (1988) [Soy . Phys. JETP 67, 1163 (1988)]. 

Y. Oreg, P. Browner, B. Simons and A. Altland, |cond-mat/9902013[ 

L. S. Levitov and A. V. Shytov, Pis'ma ZhETF 66, 200 (1997) [JETP Lett. 66, 214 (1997)]; extended version: |cond- 



mat/9607136 



F. W. J. Hekking and Yu. V. Nazarov, Phys. Rev. B 49, 6847 (1994). 

Yu. Oreg and A. M. Finkelstein, Phys. Rev. Lett., in press. 

M. V. Feigel'man and A. I. Larkin, Chem. Phys. 235, 107 (1998). 

B. L. Altshuler, A. G. Aronov and A. Yu. Zyuzin, ZhEFT 86, 709 (1984) [JETP 59, 415 (1984)]. 



39 



